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1. INTRODUCTION 


Let us consider the discrete Schrodinger equation on the half line 

-V’n+l + 2 ^n - V’n-l + V n l)) n = \lf n , U > 1, (1.1) 

where A is the spectral parameter, n is the spacial independent variable taking positive 
integer values, and the subscripts are used to denote the dependence on n. Thus, ip n 
denotes the wavefunction 4(A ,n). Note that 4o appears in (1.1) when n — 1, and hence 
the value at n — 0 of the wavefunction must be compatible with (1.1). By V n , which 
can also be written as V(n), we denote the value of the potential at n. The potential V is 
assumed to belong to class Ab, which is specified in the following definition. 

Definition 1.1 The potential V belongs to class Ab if the V n -values are real and the support 
of the potential V is confined to the discrete set n e (1, 2,. .., b} for some positive integer 
b, i.e. V n — 0 for n > b. 

As indicated in Definition 1.1, we assume that the potential in (1.1) is compactly 
supported, i.e. it vanishes for n > b + 1. Thus, the knowledge of the potential V in class 
Ab is equivalent to the knowledge of the ordered set { V\ . U 2 , • • • • 14}- We do not necessarily 
assume that all the values in this ordered set are nonzero. 

Note that (1.1) is the discrete analog of the Schrodinger equation on the half line 

—'ip"(k,x) + V(x)'ip(k,x) = k 2 if(k,x), x > 0, (1.2) 

where A := k 2 , the prime denotes the ^-derivative, the potential V is real valued, and the 
support of the potential is confined to the interval x G [0, b\. 

We are interested in the inverse problem of recovery of the potential V in (1.1) when 
the available input data set consists of the so-called transmission eigenvalues. The precise 
meaning of transmission eigenvalues is given in Section 5. The transmission eigenvalues 
are uniquely determined by the potential V and the Dirichlet boundary condition (2.2) 
specified in Section 2. Informally speaking, the transmission eigenvalues correspond to 
certain A-values at which the scattering from (1.1) with the Dirichlet boundary condition 
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(2.2) agrees with the scattering when the potential is identically zero. This inverse problem 
was recently analyzed in [22], and it was shown that the transmission eigenvalues uniquely 
determine the potential with the exception of one case [22], which we now call the “unusual 
case,” namely when the sum of the transmission eigenvalues is equal to 4(6 — 1), where 6 
is the positive integer appearing in Definition 1.1 related to the support of the potential. 

A recursive procedure was presented in [22] to determine the potential from the trans¬ 
mission eigenvalues in the “usual case,” and the nonuniqueness in the unusual case was 
illustrated via an example in [22], However, the procedure of [22] to construct the potential 
is not as direct as either of the Marchenko procedure and the Gel’fand-Levitan procedure 
we present in this paper because it requires first the construction of the right-hand side of 
(2.70) in our paper, namely the regular solution to (1.1) as a function of A, and then using 
that regular solution in (1.1) to extract the values of the potential. 

We remark that, when the potential V vanishes for n > 6+1 while V), ^ 0, as indicated 
in Theorem 5.2, the number of transmission eigenvalues (including multiplicities) is 26 — 2, 
and hence the inverse problem dealing with transmission eigenvalues is meaningful only 
when 6 > 2. 

One of our main goals in this paper is to develop a comprehensive approach to solve 
the inverse problem in the usual case in order to determine the potential V in class Ab from 
the set of corresponding transmission eigenvalues. This is done via the Marchenko method 
and also the Gel’fand-Levitan method for (1.1). Another goal we have is to elaborate on 
the unusual case and to show that the inverse problem in the unusual case may or may 
not have a unique solution. This is done by presenting various examples in Section 6 
and by showing that all of the following possibilities occur. There may be only one real¬ 
valued potential corresponding to a given set of transmission eigenvalues, there may be a 
finite number of distinct real-valued potentials corresponding to a given set of transmission 
eigenvalues, or there may be infinitely many distinct real-valued potentials corresponding 
to a given set of transmission eigenvalues. 

The literature on transmission eigenvalues and related inverse problems is expanding 
rapidly. We refer the reader to [1,3,4,7,8,12-15,19-21] and the references therein for a 
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historical account of transmission eigenvalues and some important developments in the 
field. Since we are dealing with the difference equation ( 1 . 1 ) rather than the differential 
equation (1.2), the most relevant reference for our paper is [22], The techniques used in 
our paper are analogs of the techniques from [3,4] used in the continuous case for (1.2). 

For the inverse problem for the discrete Schrodinger equation and the related Gel’fand- 
Levitan method, we refer the reader to [9], even though the discrete Schrodinger equation 
used in [9] differs from (1.1). In [9], the discrete Schrodinger equation is treated as the 
Jacobi equation 

I e -(v„+v n+ i)/ 2 ^ +i + 1 e -(v n - 1 +v n )/ 2 ^_ i = A 0 n , n > 1 . ( 1 . 3 ) 

To help the reader, in presenting the results for (1.1) we at times state the corresponding 
results for ( 1 . 2 ) by illustrating the similarities and differences between the discrete and 
continuous cases. 

Our paper is organized as follows. In Section 2 the preliminaries are presented for the 
Marchenko method to be given in Section 3, the Gel’fand-Levitan method to be given in 
Section 4, and the inverse problem with transmission eigenvalues to be given in Section 5. 
This is done by introducing the so-called Jost solution f n to (1.1), the regular solution ip n 
to (1.1), the Jost function /o corresponding to the value of the Jost solution at n = 0, and 
the scattering matrix S given in (2.45), and by presenting certain relevant properties of 
such quantities because those properties are needed in later sections. Note that, as in the 
continuous case, we call S a “matrix” even though it is scalar valued. In Section 3 the 
linear system of Marchenko equations for (1.1) with the boundary condition (2.2) is derived 
and it is shown how the potential can be recovered from the solution to the Marchenko 
system (3.12). The Marchenko system uses as input the scattering matrix and the bound- 
state information consisting of bound-state energies and the so-called Marchenko bound- 
state norming constants. In Section 3 it is also shown that, as a result of the compact- 
support property of the potential, the bound-state information is actually contained in the 
scattering matrix and hence the scattering matrix alone uniquely determines the potential. 
In Section 4 the Gel’fand-Levitan method for (1.1) with the boundary condition (2.2) 
is derived and it is shown how the potential can be recovered from the solution to the 
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linear system of Gel’fand-Levitan equations (4.17). It is also shown how the kernel of the 
Gel’fand-Levitan system can be constructed from the Jost function alone. In Section 5 the 
transmission eigenvalues related to (1.1) with the boundary condition (2.2) are introduced 
and they are shown to correspond to the zeros of the key quantity D given in (5.2) or 
equivalently to the zeros of the quantity E defined in (5.7). It is also shown how the 
Jost function / 0 can be recovered from the key quantity E in the usual case, and hence, 
having /o at hand, one can use either the Marchenko method or the Gel’fand-Levitan 
method to recover the potential. In Section 5 it is further indicated what happens when 
the unusual case occurs, which may prevent the recovery of a unique potential from the 
set of transmission eigenvalues. Finally, in Section 6 some examples are presented to 
illustrate the unique recovery of the potential from the transmission eigenvalues via the 
Marchenko method and the Gel’fand-Levitan method in the usual case. In that section 
some further examples are presented in the unusual case to illustrate that there may exist a 
unique real-valued potential for some set of transmission eigenvalues, there may be a finite 
number of distinct real-valued potentials for some other sets of transmission eigenvalues, 
or there may be infinitely many distinct real-valued potentials corresponding to some sets 
of transmission eigenvalues. 

2. PRELIMINARIES 

In this section we introduce the Jost solution f n to the discrete Schrodinger equation 
(1.1), the regular solution ip n to (1.1), the corresponding Jost function / 0 , the related 
scattering matrix S, and we present some properties of these quantities needed later on for 
the solution to the inverse problem of recovery of the potential in (1.1) from the so-called 
transmission eigenvalues. 

Using the shift operators P + and P~ and the identity operator /, which are defined 
as 

P + ^ n =y !> n +l, P-'lp n = ‘lpn-1, tyn = 
let us write (1.1) in the operator form as C^ n = A^ n , where we have defined 

C —P + + 2/ — P~ + V n . (2.1) 
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Let us associate with (1.1) the Dirichlet boundary condition 

'ipo = 0 , ( 2 . 2 ) 

which is the analog of the Dirichlet boundary condition for (1.2), namely 

ip(0) = 0. (2.3) 

Note that £ given in (2.1) with the boundary condition (2.2) is a bounded operator 
because the shift operators P + and P~ appearing in (2.1) each have norm one and the 
multiplication by V n is a bounded operator. Furthermore, £ is selfadjoint [23], which can 
be verified by establishing the equality 

oo oo 

n= 1 n =1 

where and are two square-summable sequences satisfying (2.2). Note 

that we use an asterisk to denote complex conjugation, and we recall that the square 
summability of is described by 

OO 

^ IV’nl 2 < +oo. (2.4) 

n= 1 

The bound states of £ correspond to the discrete spectrum of £, i.e. the set of those A- 
values for which there exists a corresponding square-summable solution to (1.1) satisfying 
( 2 . 2 ). 

When the potential V belongs to class Ab specified in Definition 1.1, as we will see in 
Theorem 2.5, the bound states can only occur when A G (— oo, 0) U(4, Too) and the number 
of bound states is either zero or a finite number. We will use N to denote the number 
of bound states and assume that the bound states occur at A = fij for j = 1,... ,N. The 
continuous spectrum of the discrete Sclirodinger operator £ with the Dirichlet boundary 
condition (2.2) is the interval A G [0,4] because as we will see, for each A G [0,4], we have 
a solution to (1.1) satisfying (2.2) in such a way that such a solution is either bounded in 
n without being square summable or it grows linearly in n as n —> Too. 
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In case V n = 0 in (1.1), let us use -0 n 1° denote the corresponding wavefunction. Thus, 

O 

0 n satisfies the unperturbed discrete Schrodinger equation 

o o o o 

-0n+l +20n -0n-l = A0 n , n > 1. (2.5) 

It is convenient at times to use not A but another spectral parameter related to A, 
which is usually denoted by 2 . This is done by letting 

Z := 1 - 5 + 1 vTXV T), ( 2 . 6 ) 

where the square root is used to denote the principal branch of the complex square-root 
function. Let us use T to denote the unit circle \z\ — 1, use T + for the upper portion of 
T given by z — e ie with 0 G (0, 7r), and use T+ for the closure of T + given by z — e ld 
with 6 G [0, 7r] . Under the transformation A G 2 , the point X — oo corresponds to z — 0 
in the complex 2 -plain', the point A = 0 corresponds to z — 1, and the point A = 4 to 
z = —1. The real half line A G (—oo, 0) is mapped to the real interval z G (0,1) in the 
complex ^-plane, the real interval A G (4, +oo) is mapped to the real interval 2 G (—1,0), 
and the real interval A G (0,4) is mapped to T + . Hence, the real A-axis is mapped to the 
boundary of the upper half of the unit disc in the complex 2 -plane. Thus, the bound states 
in the complex 2 -plane can only occur in the interval 2 G (—1, 0) or 2 G (0,1). From (2.6) 
it follows that 

= 1 -1 - \ vF!. (2.7) 

We see from (2.6) and (2.7) that 

A = 2-2-2 -1 . (2.8) 

Let us express 2 in terms of still another spectral parameter 6 as 

2 = e id . (2.9) 

Letting 

cos 6 = 1 -, sin 9 = - \J A(4 — A), (2.10) 

2 2 


7 



as we have already observed, as 0 takes values in the interval 6 e [0,7r] the spectral 
parameter z traces T+ and from (2.10) we see that the spectral parameter A then moves 
in the interval A G [0,4] from the left endpoint A = 0 to the right endpoint A = 4. Note 
that, from (2.8) and (2.9), we have 

A = 2 — 2 cos 6. (2.11) 

Using (2.8) it is convenient to write (1.1) as 

ifn+l + i’n-l = (z + Z- 1 + v n ) ll>n, Tl> 1 . ( 2 . 12 ) 

Thus, (2.5) can be written as 

t^n+l + V’n-l = {z + 2 - 1 )' 0 n , n > 1. (2.13) 

One can directly verify that z n and z~ n are solutions to (2.13) and that they are linearly 
independent when z ^ ±1. Hence, the general solution to (2.13) can be written as a linear 
combination of z n and z~ n when z ±1. In (2.52) we display two linearly independent 
solutions to (2.13) when z — 1. In (2.56) we list two linearly independent solutions to 
(2.13) when 2 = -1. 

There are certain relevant solutions to (2.12) and hence equivalently to (1.1). One 
such solution is the so-called regular solution ip n satisfying the initial conditions 

<p 0 = 0, ^i = l, (2.14) 

and it is the analog of the regular solution <p(k,x) to the Schrodinger equation (1.2) with 
the initial conditions 

<p(k, 0)=0, if\k,0) = 1. (2.15) 

Note that k and —k appear in the same way in (1.2) and (2.15), and hence the regular 
solution <p(k,x ) satisfies (p(—k,x ) = tp(k,x). Similarly, z and z -1 appear in the same way 
in (2.12) and (2.14), and hence tp n remains unchanged if we replace z with z~ x in <p n . 

Another relevant solution to (2.12) and hence also to (1.1) is the Jost solution f n 
satisfying the asymptotic condition 


fn — z n [l + o(l)j, 


n —> Too. 


( 2 . 16 ) 



It is the analog of the Jost solution f(k,x ) satisfying (1.2) with the asymptotics 


f(k, x) — e lkx [l + o(l)j, x —> +oo. 


(2.17) 


Recall [4,5,10,17,18] that (1.2) has also the solution g(k,x ) satisfying the asymptotics 
g(k, x) — e~ lkx [l + o(l)], x —> +oo, 

and g(k, x) is related to the Jost solution as 

g{k,x) = f{—k, x). 

In the discrete case, the corresponding analogous solution to (2.12) or equivalently to (1.1) 
is denoted by g n and it satisfies the asymptotics 


g n = £ n [l + o(l)], n +oo. 


(2.18) 


From (2.12), (2.16), and (2.18) it follows that g n is obtained from f n by replacing z by 
in f n . 

When V n = 0, let us use f n and g n to denote the corresponding Jost solution f n and 
its relative g n , respectively. From (2.13), (2.16), and (2.18) we see that 

fn = z n , °9n = z~ n , n > 1, (2.19) 

and from (2.13) and (2.19) we observe that the values of f n and g n at n = 0 are obtained 
as 

fo = 1, = I- (2.20) 


One can directly verify that the regular solution <p n satisfying (2.13) and the initial 
conditions (2.14) is given by 




Z — Z 


-1 ’ 


n > 0. 


( 2 . 21 ) 


Using (2.9) in (2.21), we can write <p n also as 


7 ^ n 


sin {nO) 
sin 9 


n > 0. 


( 2 . 22 ) 
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Equivalently, using (2.6) and (2.7) in (2.21) we can express ip n in terms of the spectral 
parameter A as 


Tn 


HhdMi) 



VA(A-4) 


i _\ n 

-VA(A-4)j 
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n > 0. 


(2.23) 


Let us now investigate the Jost solution f n to (2.12) when the potential V belongs to 
class Ab specified in Definition 1.1. Using V n = 0 for n > b + 1 in (2.12) we see that 


fn = Z n , n> b. 


(2.24) 


Let us emphasize that (2.24) holds at n = b as well. It is convenient to define 


m n := z~ n f n , (2.25) 

which is the analog of the Faddeev function m(k, x) related to (1.2) and expressed in terms 
of the Jost solution f(k,x ) as 

m(k,x ) := e~ lkx f(k, x). 

Note that from (2.24) and (2.25) we get 


m n = 1, n > b. 


(2.26) 


Using (2.24) and (2.25) in (2.12) we obtain 

m n -i — —z 2 m n+ i + (z 2 + V n z + 1) m n , n > 1, (2.27) 

and hence, we see that (2.27) allows us to obtain the values of m n and f n at n = 0. 

Proposition 2.1 Assume that the potential V belongs to class Ab specified in Definition 1.1 
and that Vi, ^ 0. Then: 

(a) For 0 < n < b — 1, the solution m n to (2.27) with the asymptotic condition (2.26) is 
a polynomial in z of degree (2b — 2n — 1) with coefficients uniquely determined by the 
ordered set {U n+ i, U n + 2 , • • •, 14 } of potential values as 

26—2n—1 

m n = ^2 K n(n+j)Z 3 , 0 < n < b - 1. (2.28) 

3=0 
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In particular, for 0 < n < b — 1 we have 


b 


H-nn 1 5 Kn(n+1) ^ ^ 

Kn(n+2) ^ ^ Vj Vl, 

(2.29) 

j=n+ 1 

n+l<j<l<b 


6-1 

K-n(2b—n— 2) Vb ^ ^ Vj, 

JV n (2b—n — l) Vfj. 

(2.30) 

j=n+1 




(b) For 0 < n < b — 1, the Jost solution f n to (2.12) satisfying (2.16) is a polynomial in z 
of degree (2b — n — 1 ) with the lowest-power term being z n and the highest-power term 
being VbZ 2b ~ n ~ 1 , and it is given by 

2b—2n—l 

fn= J2 Kn(n+j)Z n+J . ( 2 . 31 ) 

3 = 0 

(c) For 0 < n < b — 1, the solution g n to (2.12) satisfying (2.18) is given by 

2b—2n—l 

9 n= Y, Knin+j)*-^ ■ ( 2 . 32 ) 

3 = 0 

(d) The coefficients K nm are real valued, and for every nonnegative pair of integers n and 
to we have 


H-nm — 0, 

n > 2b — to, 

(2.33) 

K = 0 

n > to + 1. 

(2.34) 


PROOF: We obtain (2.28) by solving (2.27) iteratively and using (2.26). We then get (2.31) 
with the help of (2.25) and (2.28). Finally, (2.32) is obtained by using z -1 instead of z in 
(2.31). The coefficients K nm appearing in (2.33) and (2.34) are missing in the expansions 
in (2.31) and (2.32), and hence they can be chosen as zero. The real-valuedness of K nrn 
follows from the fact that f n and g n are complex conjugates of each other when z is on 
the unit circle T. | 

Note that, for 1 < n <6—1 in going from f n to f n -i, two new powers of z are gained, 
which is seen by displaying f n and f n ~i with powers in z in an ascending order as 

f n = z n + [V n+1 + • • • + V b ] z n+1 + • • • + 14 [K+r + • • • + H-r] z 2b ~ n ~ 2 + V b z 2b ~ n ~ l , 
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fn —1 — z n 1 + [V n + • • • + Vb] z n + • • • + V b [V n + • • • + Vfo_i] z 2b n 1 


+ V b z 2b ~ n . 


We remark that the coefficients K nm appearing in Proposition 2.1 will be used in Section 3 
in the Marchenko method for the inverse problem associated with (1.1). 

Recall that the Wronskian of any two solutions ip(k,x) and (f>(k,x) to (1.2) is defined 
with the help of a matrix determinant as 


b(’(k,x);<f>(k,x)] 


ip(k,x) (f>(k,x) 
1 p'(k,x) <f)'(k,x) 


(2.35) 


It is known [10,17,18] that the Wronskian in (2.35) is independent of x and hence is a 
function of the spectral parameter k alone. The analog of (2.35) for the Wronskian of any 
two solutions if n and (j) n to the discrete equation (1.1) is given by 


0n. 


V’n+l ^n+l 


Using (2.12) in (2.36) we get 


[V’nj 0n] 


V’n 




—'f’n- 1 + (z + Z 1 + V n ) ij) n —<f> n -1 + (z + Z 1 + V n ) (f) n 


V’n fn 

fn—1 fn—1 


(2.36) 


(2.37) 


From (2.37), by interchanging the two rows in the last determinant and using (2.36), we 
obtain 

[fn] fn\ = [fn— 1 j fn— l] i 


confirming that the Wronskian of any two solutions to (1.1) is independent of n and hence 
it is a function of the spectral parameter alone. Thus, the value of the Wronskian given in 
(2.36) can be evaluated at any n-value, e.g. at n = 0 or as n —> +oo. 

For z f ±1, the regular solution ip n can be expressed as a linear combination of f n 
and g n given in (2.31) and (2.32), respectively. Writing 


(fin fn T ft 9m 


(2.38) 
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where a and /? are some coefficients depending only on the spectral parameter z, we can 
express a and f3 in terms of Wronskians as 


1 ,9n'i V^n] n \fn'i n ] 

\9n‘,fn]’ [fn', 9n\ 


(2.39) 


Since the Wronskian of any two solutions to (1.1) is independent of n, we can evaluate the 
Wronskians appearing in the numerators in (2.39) at n — 0. and hence with the help of 
(2.14) and (2.36) we obtain 


L 9n, Vn\ = go, [,fn', Vn] = fo- (2-40) 

Similarly, we can evaluate the Wronskians appearing in the denominators in (2.39) as 
n —>■ +oo, and hence with the help of (2.16), (2.18), and (2.36) we get 

[fn‘,9n] = Z - 1 ~ Z. (2.41) 


Thus, (2.38) can be written as 

^Pn T iOO fn fo gn) • 

Z — Z 1 

Note that (2.40) is analogous to 


(2.42) 


[. f(~k,x);<p(k,x)\ = f(-k, 0), [f(k,x);<p(k,x)\ = f(k, 0), 


(2.41) is analogous to 


[f(k,x);f(-k,x)] = -2 ik, 


and (2.42) is analogous to 

<p(k, x) = ^ [f{~k, 0) f{k, x ) - f{k, 0) f{-k, x)]. 


(2.43) 


Let us remark that the property ip(—k,x) = ip{k,x) readily follows from (2.43). Similarly, 
from (2.42) we directly observe that replacing z with z~ x in (p n does not change the value 

of (fin- 
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It is known [10,17,18] that f(k, 0) is the so-called Jost function for the Schrodinger 
equation (1.2) with the Dirichlet boundary condition (2.3). Similarly, we see that fo corre¬ 
sponds to the Jost function for the discrete Schrodinger equation (1.1) with the Dirichlet 
boundary condition (2.2). The scattering matrix S(k) for (1.2) with (2.3) is defined as 

S(k) := (2.44) 

and similarly, the scattering matrix associated with ( 1 . 1 ) and ( 2 . 2 ) is defined as 

S := (2.45) 

Jo 

Some relevant properties of the Jost function fo are given in the following theorem. 

Theorem 2.2 Assume that the potential V belongs to class Ab specified in Definition 1.1 
and that 14 7 ^ 0. Then, the Jost function fo appearing in (2.f2) is a polynomial in z of 
degree 26—1, it is equal to mo appearing in (2.28) with n — 0, and it is given by 

26-1 

fo = E K 0j z j , (2.46) 

3= 0 

where the coefficients Koj are uniquely determined by the ordered set {Vi,..., 14} of po¬ 
tential values. In particular, we have 

b 

Koo = 1 , K (n =Y, V n K 02 = E V 3 V ^ ( 2 - 47 ) 

3= 1 i<«<i<6 

6-1 

7^0(26-2) = 7i 0 (2b-i) = V&. (2.48) 

3= 1 

Thus, fo — 1 is a “plus” function in the sense that it is analytic in z when \z\ < 1, 
continuous when \z\ < 1, and behaves like O(z) as z —> 0 in \z\ < 1. Similarly, the quantity 
go appearing in (2.f2) is given by 

26-1 

go=Y K °3 z ~\ ( 2 - 49 ) 

3 =0 

and hence go — 1 is a “minus” function in the sense that it is analytic in z when \z\ > 1, 
continuous when \z\ > 1, and 0(1/z) as z — » 00 in \z\ > 1. 
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PROOF: With the help of (2.25), the results for /o directly follow from Proposition 2.1. 
The expressions for the coefficients in (2.47) and (2.48) are obtained with the help of (2.29) 
and (2.30). The results for go are obtained from /o given in (2.46) by replacing z with 

z- 1 . I 

We observe from (2.46) and the first equality in (2.47) that the Jost function / 0 does 
not vanish at 2 = 0, and in fact the value of the Jost function at z — 0 is given by 1. It is 
convenient to display /o in powers of z in an ascending order as 

fo = 1 + [Vi + • • • + V b ] z T • • • + V b [Vi + • • • + V b -i] z 2b ~ 2 + V b z 2b -\ (2.50) 

Theorem 2.3 Assume that the potential V belongs to class A b specified in Definition 1.1. 
Then, the ordered set {Vi, V 2 ,..., V b } of potential values is uniquely determined by the 
ordered set {Kqi, K 12 , K 23 , ■ ■ ■, K(b-i)b}, where the quantities K nm are the coefficients 
appearing in (2.31). 

PROOF: From the second identity in (2.29) we obtain 

^(n—l)n -^n(n+l)> 1 5; ^ b, (2.51) 

where we also use K b r b+ 1 ) = 0, which follows from (2.33). Hence, the proof is complete. | 

The following theorem describes the behavior of the regular solution <p n as n Too 
when A = 0 and when A = 4. It shows that (1.1) with the Dirichlet boundary condition 
(2.2) cannot have a bound state at A = 0 or A = 4. It also shows that (1.1) and (1.2) have 
some similarities for the generic case and also for the exceptional case. 

Theorem 2.4 Assume that the potential V belongs to class Ab specified in Definition 1.1. 
Let A and z be the spectral parameters appearing in (1.1) and (2.12), respectively, and let 
ip n and f n be the corresponding regular solution and the Jost solution to (1.1) appearing in 
(2.14) and (2.16), respectively. Let f n and p n be the Jost and regular solutions appearing 
in (2.19) and (2.21), respectively, corresponding to V n = 0. Then: 

(a) At A = 0, and equivalently at z = 1, the regular solution ip n either grows linearly 
as n —»• Too, which corresponds to the “generic case, ” or it remains bounded, which 
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corresponds to the “exceptional case. ” Hence, A — 0 never corresponds to a bound 
state for (1.1) with the Dirichlet boundary condition (2.2). 

(b) At A = 4, and equivalently at z = —1, the regular solution ip n generically grows 
linearly as n —>■ +oo or in the exceptional case it remains bounded. Hence, A = 4 
never corresponds to a bound state for (1.1) with the Dirichlet boundary condition 
(2.2). 

(c) If fo has a zero at z = 1, then it must be a simple zero. If fo has a zero at z — 1, then 
go also has a simple zero at z — 1. 

(d) If fo has a zero at z — —1, then it must be a simple zero. If fo has a zero at z = — 1, 
then go also has a simple zero at z = —1. 

(e) The Jost function fo and the quantity go cannot vanish simultaneously at any z-value 
with the possible exception of z — ±1. 

(f) The exceptional case occurs when fo = Oatz — l or fo — 0 at z — — 1. We may have 
fo vanishing at z — 1 but not at z = —1, we may have fo vanishing at z — — 1 but not 
at z = 1, and we may have fo vanishing at both z = 1 and z = — 1. 

PROOF: For A = 0, we can use the method of variation of parameters [6] and express any 

solution to (1.1) in terms of two linearly independent solutions to (2.5). In particular, for 

A = 0 we can use the regular solution ip n and the Jost solution f n that are given by 

<Pn = n, f n — 1. n > 0. (2.52) 


The method of variation of parameters yields the following expression for the regular 
solution to (1.1): 


Tn ^ 


Thus, for n > b, we obtain 


n— 1 


i - E v > v, 

3 = 1 


n— 1 


+ n > 2 - 
3 = 1 


(2.53) 


(p n n 


i-E y ^ 

j=. l 


+ ^ ) V :i Tj , n>b+ 1. 
3 = i 


(2.54) 
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From (2.54) we conclude that tp n at A = 0 grows linearly in n as n —> +oo unless the 
exceptional case occurs with 

b 

£+^ = 1. (2.55) 

3 =1 

If (2.55) holds, the second summation term in (2.54) must be nonzero because otherwise 
we would have ip n — 0 for A = 0 and n > 6+1, causing p n = 0 for A = 0 and for all 
n > 0, contradicting the second equality in (2.14). Thus, in either case, i.e. whether (2.55) 
holds or not, from (2.54) we conclude that ip n at A = 0 cannot be square summable, and 
hence there cannot be a bound state at A = 0 for (1.1) with the boundary condition (2.2). 
For A = 4 or equivalently z = —1, the proof is obtained similarly, in which case (2.52) is 
replaced with 

£n = (- l) n “V fn = (~ l) n , n>0, (2.56) 


(2.53) is replaced with 


(p r 


(- 1 ) 


n —1 


n 


71—1 


i-£(-i yvin 


3 = 1 


+ (-l) ? 


n— 1 

■'EM 1 

3 = 1 


jVjVsi 


n > 2, 


(2.57) 


(2.54) is replaced with 


— ( 1 ) 


n — l rt 


n 


i - 


£(-!)+¥5- 


3 = 1 


+ (-1)”- 1 J](-1)P Kj n>6+l, (2.58) 


3 = 1 


and (2.55) is replaced with 

b 

J](-i)W iW = i. 

t=i 

Thus, the proofs of (a) and (b) are complete. Let us now prove (c), (d), and (e). Since 
go is obtained from /o by replacing z by z _1 , it follows that go vanishes at z — 1 if fo 
vanishes there. Similarly, go vanishes at z = — 1 if fo vanishes there. However, fo and go 
cannot vanish simultaneously at any other z -value because otherwise, as a result of (2.42), 
the regular solution ip n would be zero at that 2 -value for all n > 1, contradicting p\ — 1 
stated in (2.14). Thus, in order to complete the proofs of (c), (d), and (e), we only need 
to show that a possible zero of fo at 2 = 1 must be a simple zero and that a possible 
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zero of /o at z — —1 must be a simple zero. By Proposition 2.1(b) we know that f n is a 
polynomial in z and this is true even when 14 = 0. Hence, we can expand f n around z — 1 
as 

fn(z) = /n(l) + (z - 1) /n(l) + O ((z ~ l) 2 ) , Z 1 ill C, (2.59) 

where f n (z) denotes the value of f n at z and f n ( 1) denotes the derivative of f n with respect 
to z evaluated at z = 1. Since g n is obtained by replacing z with 2 _1 in f n , from (2.59) we 
get 

-l) /„(!) + O 

Using 

- — 1 = —(z — 1) + O ((z — l) 2 ) , z —> 1 in C, 

z v ' 

we can write (2.60) as 

9n(z) = /n(l) - (z - 1) /n(l) + O ((z - l) 2 ) , Z 1 ill C. (2.61) 


^n(^) = fn{ 1) + 


- 1 


z ->• 1 in C. 


(2.60) 


Using (2.36), (2.59), and (2.61), we write the left-hand side of (2.41) as 

fn( 1) fn( 1) 

/n +1 (l) /n+l(l) 


[fn(z); g n (z)\ = 2(2-1) 
Since we have 


+ 0((^-l) 2 ), z^linC. (2.62) 


= 2{z - 1) + O ((z - l) 2 ) , 2 —^ 1 in C, 

with the help of (2.41) and (2.63) we obtain 


(2.63) 


[fn(z)\g n {z)\ = -2(z - 1) + O ((z - l) 2 ) , 2 ->• 1 in C. (2.64) 


Comparing (2.62) and (2.64) we see that 

fn( 1) 


fn( 1) 
fn+1 (1) /n+l(l) 


= 1 . 


(2.65) 


From (2.65) we conclude that f n ( 1) and f n { 1) cannot vanish simultaneously. Thus, if 
f n { 1) = 0, then we must have / n ( 1) ^ 0. Therefore, 2 = 1 can only be a simple zero of 
fn(z)- Proceeding similarly, we can prove that if f n has a zero at 2 = —1, then such a zero 
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must be a simple zero. Thus, the proof of (c)-(e) is complete. Finally, let us prove (f). 
For this, let us first show that ip n at z = 1 remains bounded as n —* +oo if and only if 
/o(l) = 0. Using (2.59) and (2.61) with n — 0 as well as using (2.12), (2.13) and (2.62), 
we expand the right-hand side of (2.42) around z — 1 as n —* +oo, and we obtain 

Tn = n/ 0 (l) -/o(l) + 0(z- 1), z-¥ 1 in C, (2.66) 

where the terms represented with 0(z — 1) remain bounded as n —> +oo. From (2.66) we 
conclude that ip n at z — 1 remains bounded if and only if / 0 (1) = 0, proving that the 
exceptional case with z — 1 occurs if and only if /o(l) = 0. A similar argument proves 
that the exceptional case with z — —1 occurs if and only if /o(—1) = 0. We can show with 
some explicit examples that /o can vanish at z — 1 but not at z = —1, and vice versa. 
For example, in Example 6.1 in Section 6, from (6.1) we see that we have /o(l) = 0 and 
/o( — 1) = 2 if we choose V\ — —1. In that example, we get /o(—1) = 0 and /o(l) = 2 if 
we choose V\ = 1. In Example 6.3, choosing V\ — —y/2 and V 2 — l/\/2, from the second 
equality in (6.12) we see that /o(l) = 0 and fo(— 1) = 0. | 

We remark that the summation terms appearing in (2.57) and (2.58) are the discrete 
analogs of the integrals appearing in (5.5) and (5.6) of [2], In the next theorem, we 
summarize the facts relevant to the bound states of (1.1) with the boundary condition 
(2.2). Recall that the bound states correspond to the A-values at which (1.1) has a square- 
summable solution satisfying the boundary condition (2.2). 

Theorem 2.5 Assume that the potential V belongs to class Ab specified in Definition 1.1. 
Let A and z be the spectral parameters appearing in (1.1) and (2.12), respectively, and let 
f n , ip n , and fo be the corresponding Jost solution appearing in (2.16), the regular solution 
appearing in (2-42), and the Jost solution appearing in (2-46). Then: 

(a) A bound state can only occur when A G (— 00 , 0) or A G (4, + 00 ). Equivalently, a bound 
state can only occur when z G (—1, 0) or z G (0,1). 

(b) The regular solution evaluated at a bound state, i.e. the value of ip n evaluated at 
a X-value corresponding to a bound state must be real. Similarly, the Jost solution 
evaluated at a bound state must be real. 
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(c) At a bound state the Jost function fo has a simple zero in A and in z. At a bound state 
the value of the Jost solution at n — 1 cannot vanish, i.e. fi^Oata bound state. 

(d) The number of bound states N must be finite, and we have 0 < N < 2b — 1, where b 
is the positive number related to the support ofV and appearing in Definition 1.1. In 
particular, we have N — 0 when V n = 0. 

(e) At a bound state the quantity go appearing in (2.f9) cannot vanish. Thus, the scat¬ 
tering matrix appearing in (2.f5) has a simple pole at each bound state. 

PROOF: By its definition, at a bound state, (1.1) must have a square-summable solution 
and that solution must satisfy (2.2). Since the regular solution p> n satisfies (2.2), any 
bound-state solution must be linearly dependent on the regular solution. Since the cor¬ 
responding operator £ is selfadjoint, a bound state can only occur when A is real valued. 
Recall that, as seen from (2.6), the real A-values in the interval (0,4) correspond to the 
z- values on T + , the upper portion of the unit circle. On the other hand, the z- value cor¬ 
responding to a bound state cannot occur on T + , because, as a consequence of (2.16) and 
(2.18), neither of the two linearly independent solutions f n and g n to (1.1) vanish at those 
^-values as n —» +oo. Furthermore, from Theorem 2.4 we know that a bound state cannot 
occur at z — 1 or at z = — 1. We also know from (2.6) that z — 0 corresponds to A = oo and 
hence a bound state cannot occur at z = 0. Thus, we have proved (a). Let us now prove 
(b). Since the bound states can only occur at some real z- values, from Proposition 2.1(d) 
we conclude that f n given in (2.31) and g n given in (2.32) are both real valued at a bound 
state for any n > 0. As a result, from (2.42) we conclude that the regular solution ip n also 
takes real values at a bound state. Thus, we have proved (b). Let us now turn to the proof 
of (c). At a bound state, we conclude that the coefficient of g n on the right-hand side of 
(2.42) must be zero because otherwise (2.18) implies that ip n cannot be square summable 
at the bound state occurring at some 2 - value lying in the interval z G (—1, 0) or z G (0, 1). 
Thus, at a bound state the Jost function / 0 must be zero. At a bound state we cannot 
have fi — 0 because then the only solution to (1.1) with /o = 0 and fi — 0 would have to 
be f n = 0, which is not compatible with (2.16). Let us now prove that the zero of fo at a 
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bound state must be simple. From (1.1) we obtain 


fn+ 1 + fn-l — (2 — A + V n ) fn , 


dfn+1 df n —1 . \dfn 

~ir + ^r~~ fn + {2 ~ x+Vn h a ’ 


(2.67) 


where the second line is obtained by taking the A-derivative of the first line. Multiplying 
the first line in (2.67) by df n /d\ and the second line by f n and taking the difference of the 
resulting equations we get 


r dfn dfn- |-i clf n 

J Tl~\- 1 j x J 71 ' J Tl 1 7 A 

aX aX aX 


~fr 


dfn-1 

dX 


= /, 


( 2 . 68 ) 


Let us evaluate (2.68) at a A-value corresponding at a bound state and take the summa¬ 
tion of both sides starting with n — 1. Many cancellations on the left-hand side in the 
summation yield 


dfi , dfp __ 2 

dX 11 r1\ ~ 2^ J ri- 


dfo 

dX 


(2.69) 


n= 1 


We already know that at a bound state /o = 0. f\ =7 0, and the / n -values are real. Thus, 
(2.69) implies that at a bound state we have 


dfo _ _ J_ \ ^ „2 

dX ~ hk n ' 

with the right-hand-side being nonzero. Thus, the zero of /o as a function of A at the 
bound state must be simple. From (2.6) we see that dX/dz at a bound state is strictly 
positive and hence a simple zero of /o in A corresponds to a simple zero of /o in z. Thus, 
the proof of (c) is complete. By Theorem 2.2 we know that /o is a polynomial in z of 
degree 2b — 1. Since the bound states correspond to the zeros of /o, we conclude that the 
number of bound states cannot exceed 26—1. Furthermore, from the first equality in (2.20) 

O 

we know that the Jost function f n corresponding to the zero potential does not have any 
zeros and hence the number of bound states corresponding to the zero potential is zero. 
Thus, we have proved (d). At a bound state z -value we already know that /o vanishes and 
that 2 -value cannot be equal to 1 or —1. Thus, from (2.42) we see that the vanishing of go 
at a bound state would imply ip n = 0, contradicting the second condition in (2.14). Then, 
from (2.45) and the simplicity of the zeros of /o at the bound states, we conclude that the 
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scattering matrix S must have a simple pole at each bound state. Thus, the proof of (e) 
is complete. | 

In (2.42) the regular solution tp n to (1.1) with the initial conditions (2.14) is expressed 
in terms of the spectral parameter z. It is possible to express ip n in terms of the spectral 
parameter A appearing in (1.1). That result is stated in the next theorem and will be 
useful in the formulation of the Gel’fand-Levitan procedure for (1.1). 

Theorem 2.6 Assume that the potential V belongs to class Ab specified in Definition 1.1. 
Then, for n > 1 the regular solution ip n to (1.1) with the initial conditions (2.14) is a 
polynomial in A of degree (n — 1) and is given by 

71— 1 

<p n = J2 B njX, (2.70) 

3=0 

where the coefficients B n j are real valued and are uniquely determined by the ordered set 
{V\, V 2 , ..., 14.-i} of potential values. In particular, we have 

71—1 

B n[n - i ) = (-l) n " 1 , B n{n _ 2) = (-l)- 2 2(n~l) + J2 V i ’ ( 2 - 71 ) 

l=i 

71—1 

B n{n - 3 ) = (-l) n “ 3 (n-2)(2n-3) + 2(n-2)^^-+ ^ V 3 V k . (2.72) 

j =1 l<k<j<n —1 

PROOF: As seen from (1.1), <p n satisfies 

Tn +1 + An -1 = (2 — A + 14) (p n , n > 1. (2.73) 

Solving (2.73) iteratively and using the initial values (fo — 0 and pi — 1 as stated in (2.14), 
we get (2.70) with the coefficients given in (2.71) and (2.72). Since the initial values are 
real valued and the coefficients in (2.73) for p n -\ and ip n are real valued for real A-values, 
the coefficients B n j appearing in (2.70) are all real valued. | 

With the help of Theorem 2.6 we see that, for n > 1, the regular solution ip n appearing 
in (2.23) is a polynomial in A of degree n — 1. The same result can also be obtained directly 
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from (2.23) by using the expansions 


1 - o ± o \/-M A - 4) ) = 


3=0 


A 


n-j 


^±- a/A(A - 4)^ , 


n > 1, 
(2.74) 

with ('!) denoting the binomial coefficient n\/(j\ (n — j)!). With the help of (2.74), from 
(2.23) we obtain 


L(n—l)/2j 

Fn = 

s=0 


n 

2s+ 1 


1 


A 


n— 1 — 2s 


A S | j — 1 

4 


n > 1, 


(2.75) 


with |_^J denoting the floor function, i.e. the greatest integer less than or equal to x. 

Using either (2.75) or Theorem 2.6, we directly obtain the following corollary. 

Corollary 2.7 For n > 1, the regular solution to (2.5) with the initial conditions (2.14), 
namely ip n appearing in (2.23), is a polynomial in A of degree (n — 1) and is given by 


n — 1 q 

ip n ^ ^ 11 nj X' . 
3=0 


(2.76) 


where the coefficients B n j are real valued and obtained from the expansion in (2.75). In 
particular, we have 


B«(n- 1 ) - (-l) n “\ B n{n _ 2 ) = 2( —l) n-2 (n - 1), 
B n (n- 3 ) — (—l) n-3 (n — 2)(2n — 3). 


(2.77) 

(2.78) 


Let us now consider expressing the regular solution cp n to (1.1) in terms of the elements 
in the set {pj}f =1 . This problem arises in the formulation of the Gel’fand-Levitan method, 
and it is also closely related to the theory of orthogonal polynomials [9]. 

Theorem 2.8 Assume that the potential V belongs to class Ab specified in Definition 1.1. 
Let (f n be the regular solution to (2.5) satisfying the initial conditions (2.14). Then, for 
n> 1 the regular solution (p n to (1.1) with the initial conditions ( 2 . 14 ) can be written as 
a linear combination of the elements in the set {<pj}” =1 , and we have 

n— 1 

Fn — Fn T ^ ^ A n j ipj , n 1, (2.79) 

3 = 1 
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with the coefficients A n j for j = — 1 uniquely determined by the ordered set 

{Vi,, V n -\}. In particular, we have 

n —1 

A n (n-1) = V i' n - 2 ’ ( 2 - 80 ) 

3= 1 


with the convention that 


A-nO 0, ^4( n +l)(n+l) 1) U ^ 0, 

A nm = 0, 0 < n < m. 


(2.81) 


PROOF: From Theorem 2.6 we see that cp n is a polynomial in A of degree n — 1. From 
Corollary 2.7 we know that, for j > 1, the quantity (fj is a polynomial in A of degree j — 1. 
Thus, we can use the set {ipj}f =1 as a basis for polynomials in A of degree n — 1, and we 
express the polynomial ip n as 

n— 1 

ip n = A nn ip n T ^ ^ -A-nj ^Pjt n 1. (2.82) 

i=i 

With the help of (2.70) and (2.76), we can compare the coefficients of A n_1 on both sides 

o 

of (2.82). From the value of B n r n _ p in (2.71) and the value of £> n ( n -i) in (2.77) we know 
that those coefficients are both equal to (—l) n_1 , and hence we must have A nn = 1 in 
(2.82). Let us remark that we can use A n0 for n > 1 as the coefficient of <po 5 by recalling 
that <po — 0 as a result of the Dirichlet boundary condition given in (2.14). Thus, the 
choices in (2.81) are appropriate. Using (2.70), (2.76), and an analog of (2.76) but with 
( n — 1) instead of n there, we compare the coefficients of X n ~ 2 on both sides of (2.82). This 
leads us to the equation 

O O 

B n (n— 2) A n nB n (n— 2) T A n ^ n _ i)_B( n _i)( n _ 2 ) , Tl ^ 2. (2.83) 

Inserting in (2.83) the value A nn — 1 and the value of B n ( n _ 2 ) from (2.71) and the values 

o o 

of R n(n _ 2 ) and R (n -i)( n _ 2 ) from (2.77), we get (2.80). | 

We remark that the coefficients A nm will appear in Section 4 as the unknown quantities 
in the Gerfand-Levitan system for (1.1) with the Dirichlet boundary condition (2.2). From 
(2.80) and the first equality in (2.81) we obtain the following corollary. 
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Corollary 2.9 Let V be a potential in class Ab specified in Definition 1.1. Then, V is 
uniquely determined by the ordered set {A21, A32, • • •, ^4(b+i)b} of the corresponding coeffi¬ 
cients appearing in (2.79) via 

Vn ^-(n+l)n -^-n(n—1)> ^ ^ (2.84) 

with the understanding that Aiq — 0. 

3. THE MARCHENKO EQUATION 

In this section we outline the Marchenko method to recover the potential V from 
the corresponding scattering data for (1.1) with the Dirichlet boundary condition (2.2). 
We also show that the bound-state information is contained in the scattering matrix S 
appearing in (2.45) as a result of the compact support of V. Thus, the corresponding 
scattering matrix alone determines the potential V in class Ab uniquely without any need 
to specify the bound-state information separately. 

The Marchenko equation for (1.2) is given by [5,10,17,18] 

/»oo 

K(x, y) + M(x + y) + / dt K(x,t) M(t + y), y > x, (3-1) 

J X 

in which the unknown function K(x,y ) is related to the Jost solution f(k,x ) appearing in 
(2.17) as 

i r°° 

K(x,y):= — f dk[f(k,x)-e ik *]e- ik \ 

or equivalently as 

i r°° 

K(x,y) := — / dk{f(-k,x)-e~ ik * 

27T J_ oq 

Note that K(x, y) — 0 for x > y as a result of the facts [5,10,17,18] that if the potential V 
in (1.2) is real valued and integrable on x G (0, b ) and vanishes for x > b. then for each fixed 
x > 0 the Jost solution f(k,x ) is analytic in k in the open upper-half complex plane C + , 
continuous in k in the closed upper-half complex plane C+, and e~ lkx f(k, x) = 1 + 0(l/k) 
as k —» oo in C+. The kernel of the Marchenko equation is expressed in terms of the 
scattering matrix S(k) given in (2.44) and the bound-state information as 

i /»oo N 

M(y):=— dpi-SWJe^ + ^cP-*-", (3,2) 
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where the values k — iK a for s = 1,..., N with distinct positive k s correspond to the zeros 
of /(/c, 0) and the constants c s denote the Marchenko bound-state norming constants given 
by 

1 

c s := , ^=- 

dx f(in s ,x) 2 

Because the potential V is assumed to be supported within the finite interval x G [0, 6], it 
turns out that the value of the norming constant c s is uniquely determined by the scattering 
matrix S(k) as [3,4] 

c s — \Ji Res[5'(/c), in s \, (3-3) 

with Res[<S(fc),i«: s ] denoting the residue of S(k) at the pole k = in a . We note that the 
poles of S(k) in C + are all simple and that the number of bound states denoted by N is 
either 0 or a positive integer. 

One can derive (3.1) from (2.43) as follows. From (2.43) we get 

f(k,0)f(-k,x)-f(^k,0)f(k,x) = -2ik(p(k,x), k e R, (3.4) 

where R denotes the real axis. Dividing (3.4) by f(k, 0) and using (2.44) we obtain 

2 jk 

f(—k, x) - S(k) f(k, x) = - — <p{k,x), ke R. (3.5) 

f{k, 0) 

We can write (3.5) as 

[f(-k, x) - e~ ikx ] + [1 - S{k)] e ikx + [1 - S(k)\ [f{k , x) - e ikx ] = H{k , x), ke R, (3.6) 

where we have defined 

o jk 

H{k, x) := f(k, x) - e~ lkx - <p(k, x). (3.7) 

j(M) 

With the help of the Fourier transform of (3.6) with (27r) _1 dk e lky for y > x, we 
obtain (3.1). When the potential V in (1.2) is real valued and compactly supported, the 
bound-state information can be obtained [3,4] from the scattering matrix S given in (2.44) 
because in that case S(k) has a meromorphic extension from k e R to k e C + with simple 
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poles at k = in s for s — 1,..., N and the residues of S(k) at k — in s are related to the 
the bound-state norming constants as in (3.3). Note that the contribution of the Fourier 
transform of (3.7) from each bound state at k ~ in s is evaluated with the help of (3.4) and 
the fact that f(in s , 0) = 0, and that contribution is given by 


1 

27T 



2 ik 

7(M) 


<p(k,x) e iky = 


f(-iK a , 0) 

/(*««, 0 ) 


f(in a ,x)e KsV , 


where an overdot denotes the derivative with respect to k. From (2.44) we see that 


Kes[S(k), in s \ — 


f(-in s , 0) 


/(*««, 0 ) 

and hence c s appearing in (3.2) is given by (3.3). 

The analog of (3.1) for (1.1) is derived in a similar way. We can write (2.42) as 


(3.8) 


So 9n ~ gtlfn ~ -(« “ « ') Vn, 2 6 T, 


(3.9) 


where we recall that T denotes the unit circle \z\ — 1 in the complex z- plane. Dividing 
(3.9) by the Jost function /o and using (2.45) we get 

9n - Sfn = - ^Pm ZE T. (3.10) 

Jo 

We can write (3.10) as 


(9n - z~ n ) + (l-S)z n + (l-S)(f n -z n ) = H n , ze T, (3.11) 


where we have defined 

H n '. = f n z~ n “ *p n - 

Jo 

Let us take the Fourier transform of (3.11) with (27ri) _1 f clz z m ~ 1 for m > n + 1, where 
the integral is along the unit circle T in the counterclockwise direction traversed once. 
This yields the linear system of Marchenko equations 

2b—n— 1 

-K-nm E K n jMj +m = 0, 0 < n < m, (3-12) 

j=n +1 
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where we have defined 


K nm ■= dz (g n - z n )z m \ (3.13) 

i r N 

M n := — & dz (1 - S ) z n ~ x + y "c 2 s z n s . (3.14) 

2m J ' 

J=S 

In (3.14), the values z = for s = 1,..., N correspond to the zeros of /o inside the unit 
circle T and the Marchenko bound-state norming constants c s are defined as 


1 



(3.15) 


where we have used the fact that the value of f n at a bound state is real as asserted in 
Theorem 2.5(b), and hence we have replaced the absolute square in the denominator in 
(3.15) with the ordinary square. Since the potential V is compactly supported, it turns 
out that the Marchenko norming constant c s defined in (3.15) can be evaluated from the 
scattering matrix S appearing in (2.45) via the residue of S/z at the bound-state pole 


z — z s . We have 


/ Res 

~S ' 

— ,z s 

/ 

z 


s = l,...,N. 


(3.16) 


We recall that each bound state z — z s occurs when either z s e ( — 1, 0) or z s G (0,1). Lot¬ 
us remark that (3.13) can also be written as 


K 


nm 


l 

2ni 


dz(f n -z n )z~ m -\ 


(3.17) 


and that the contribution to the Fourier transform from the right-hand side of (3.11) at 
each zero z = z s of the Jost function / 0 is evaluated by using 


1 

2ni 


dz H„ z 


771—1 


1 


z — z 


-l 


=-® dz - 

2m J 

zfo 

_ 90 r 

z Jn 

JO 



VnZ 


(3.18) 


where we have used (3.9) and the fact that /o vanishes at z — z s . The overdot in (3.18) 
indicates the derivative with respect to z. The value of c s appearing in (3.14) is given by 
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(3.16) by noting that from (2.45) and (3.18) we have 


Res 


'S 

Z 



which is the analog of (3.8) in the discrete case. 


(3.19) 


Let us note that the quantities K nm given in (3.13) and (3.17) coincide with those 
appearing in (2.28)-(2.32) and (2.46)-(2.49). The following theorem states that the ordered 
set {14, ...,V b } of potential values is uniquely determined by the scattering matrix S 
appearing in (2.45) as well as by the Jost function / 0 appearing in (2.46). 


Theorem 3.1 Assume that the potential V belongs to class Ab specified in Definition 1.1. 
Then, the Jost function fo appearing in (2.46) uniquely determines the potential. Similarly, 
the scattering matrix S given in (2.45) uniquely determines the potential. The recovery 
of the potential can be accomplished by solving the Marchenko system (3.12) for K nm and 
then obtaining the potential V via (2.51). 

PROOF: Given the Jost function fo as a function of z, we also have go given in (2.49) 
because go is obtained by replacing z with z _1 in fo- Thus, as seen from (2.45) we have 
at hand the scattering matrix S' as a function of the spectral parameter z. On the other 
hand, if we are given the scattering matrix S, using (2.45) and (3.16) we then construct 
the Marchenko kernel M n given in (3.14). Next, we solve the Marchenko system (3.12) 
and uniquely obtain K nm for m > (n + 1) for n — 0,1,..., b. Finally, we use (2.51) to 
obtain the ordered set {14,..., 14} of potential values. | 


4. SPECTRAL FUNCTIONS AND THE GEL’FAND-LEVITAN METHOD 

Let us use p to denote the spectral function [9] associated with the unperturbed 
discrete Schrodinger operator corresponding to (2.5) with the Dirichlet boundary condition 
(2.2). Similarly, let us use p to denote the spectral function associated with the perturbed 
discrete Schrodinger operator (1.1) with the Dirichlet boundary condition (2.2). Our goal 
in this section is to derive the corresponding Gel’fand-Levitan system and establish the 
recovery of the potential V in class Ab from the spectral function p. Recall that as the 
spectral parameter we can use A appearing in (1.1), or use z appearing in (2.6), or even 
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use 9 appearing in (2.9). Thus, p and p can be viewed as functions of A or functions of z 
or functions of 9. We will refer to dp and dp as the spectral measures associated with (2.5) 
and (1.1), respectively, with the Dirichlet boundary condition (2.2). 


As seen from (1.1) and (1.3), we observe that the spectral parameter used in [9] is not 
the same as the spectral parameter A used in (1.1) but the two are closely related to each 
other. Furthermore, the discrete Schrodinger equation used in [9] is not the same as (1.1) 
we use. Nevertheless, proceeding as in [9], we obtain our spectral measure dp in terms of 
A as 


'0, A < 0, 

dp = -L a/A(4 - A) dX, 
Zir 


0 < A < 4, 


0, 


A > 4. 


(4.1) 


As a function of 6 , the nonzero portion of dp corresponding to 0 < A < 4 can be expressed 
as 

dp = — sin 2 9 dO , 0 < 9 < it. (4.2) 

7T 

As a function of 2 , the nonzero portion of dp corresponding to 0 < A < 4 is given by 

4=~(z-z- ‘) 2 -, Z6T+, (4,3) 

2m z 


where we recall that T+ denotes the closure of the upper portion of the unit circle T. We 
already know from Theorem 2.5(a) that (2.5) with the Dirichlet boundary condition (2.2) 
does not have any bound states and hence the expressions for dp do not contain any terms 
related to bound states. 

The idea [9,11,16,18] behind the spectral function p is the following. The regular 
solution tp n appearing in (2.21)-(2.23) as a function of n forms the sequence {<Pn} < £L\, 
where each term is a function of the spectral parameter A. In fact, by Theorem 2.6 we 
know that, for n > 1, the quantity <p n is a polynomial in A of degree n — 1. The terms in 
the sequence {< °p n }^ =1 form an orthonormal set with respect to the spectral measure dp, 
i.e. we have 

J dp (£j £j) = 5 jt , (4.4) 
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where Sji is the Kronecker delta and we omit the integration limits for notational simplicity 
with the understanding that the integral is over A G R. 


The orthonormality condition given in (4.4) can be stated in A, in or in 9. For 
example, in terms of 9 , as seen from (2.22) and (4.2) we can write the left-hand side of 
(4.4) as 


dp {^p 3 ipi'j = — J d6 sin 2 6 


sin (j 0) 

sind 


sin(Z 9) 
sind 


Since we already know that 


7r 


dO sin(j 9) sin(Z 9) = - 5 jt , 


we see that (4.4) holds when it is expressed in 9. 

Let us now consider the main idea behind the spectral function p associated with 
(1.1) with the Dirichlet boundary condition (2.2). The regular solution <p n to (1.1) with 
the initial conditions (2.14) as a function of n forms the sequence where each 

term in the sequence is a function of the spectral parameter A. In fact, by Theorem 2.6 we 
know that for n > 1 the quantity <p n is a polynomial in A of degree n — 1. The terms in the 
sequence {(p n }^L 1 are orthonormal with respect to the spectral measure dp , i.e. we have 


dp (c Pj <pi) Sji, 


(4.5) 


where again for simplicity we suppress the integration limits in our notation. 

The spectral function p contains all the relevant information about the potential V 
and the boundary condition (2.2). In the presence of the potential V, assuming that the 
bound states occur at A — p s for s — 1,..., N we have 

( N 


dp = { 


Cg S (A — p s ) dX, A < 0 or A > 4, 
dX 


S = 1 

1 


(4.6) 


2 Ir ' /A(4 A) |/o| 2 ’ 


0 < A < 4, 


where /o is the Jost function appearing in (2.46). Notice that the continuous part of dp in 
(4.6) is obtained from dp given in (4.1) via a division by |/o| 2 - 
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Recall that, as stated in Theorem 2.5(c), the bound states correspond to the zeros of 
the Jost function /q. The constant C s appearing in (4.6) is the Gel’fand-Levitan bound- 
state norming constant at the bound state X — p s , i.e. 



where we have used the fact that the value of <p n at a bound state is real, as stated 
Theorem 2.5(b). 

The contribution to the spectral function p from each eigenvalue is given by the 
square of the norming constant C s appearing in (4.7). We can obtain the spectral density 
given in (4.6), and in particular its continuous part, as the limiting case p —> Too from 
the spectral density p p corresponding to the discrete Schrodinger equation on the finite 
lattice n G {1,2,...,_p} with the boundary condition (2.2) at n = 0 and an additional 
boundary condition at n = p + 1, e.g. the Dirichlet condition ip p+ 1 = 0. The discrete 
Schrodinger operator corresponding to such a problem on the finite lattice has only the 
discrete spectrum and we can explicitly evaluate the corresponding spectral function p p . 
We then get the continuous part of the expression in (4.6) by taking the limit of p p as 
p —>• Too. 

By comparing the Marchenko norming constant c s defined in (3.15) with the Gel’fand- 
Levitan norming constant C s defined in (4.7), we observe that they differ from each other 
in the sense that the Marchenko norming constant c s is obtained by normalizing the Jost 
solution f n at the bound state A = p s whereas the Gel’fand-Levitan norming constant is 
obtained by normalizing the regular solution <p n at the same bound state. On the other 
hand, we already know that the regular solution p n and the Jost solution f n are linearly 
dependent at a bound state. Thus, evaluating (2.42) at the bound state A = p s and using 
the fact that /o vanishes at A = /j. s we obtain 

Since the quantity inside the parentheses in (4.8) is independent of n, by squaring both 
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sides of (4.8) and then using a summation over n, we obtain 



Using (3.15) and (4.7) in (4.9) we conclude that the Marchenko and Gel’fand-Levitan 
bound-state norming constants are related to each other as 


— c„ (4.io) 

^ Z A=/x s 

If we know |/ 0 |, because go is obtained from / 0 by replacing z with z _1 , we can evaluate 
the factor in (4.10) relating the Marchenko norming constant c s and the Gel’fand-Levitan 
norming constant C s . We remind the reader that z and z~ x are related to A as in (2.6) 
and (2.7), respectively. 



As already mentioned, the continuous spectrum of the discrete Schrodinger operator 
associated with (1.1) and (2.2) is A G [0,4] and the discrete spectrum consists of at most 
a finite number points A = p s where each /j s is located either in the interval A G (—oo, 0) 
or A G (4, +oo). Viewed as a function of the spectral parameter 0. the continuous part of 
the spectral measure dp is given by 

2 dQ 

dp = — sin 2 6 -rr~\o , 0 < 9 < 7r, (4.11) 

7T I/or 

which is obtained from (4.6) with the help of the transformation 9 i —> A given in (2.11). 
Similarly, we obtain the same quantity as a function of the spectral parameter z by using 
(2.8) in (4.6) and obtain 

1 d T' _ 

Note that, in (4.6), (4.11), and (4.12), instead of |/o| 2 , we can equally use the product foQo 
because |/o| 2 = foPo when z is on the upper unit semicircle T + because g 0 is obtained 
from /o by replacing z by z~ l . which is equivalent to taking the complex conjugate of /o 
when z is on T + . 


We derive the linear system of Gel’fand-Levitan equations as follows. Consider the 
regular solution <p n to (1.1) with the initial conditions (2.14). The set }y=i forms an 
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orthonormal basis for polynomials in A of degree n — 1 with respect to the spectral measure 
dp given in (4.6). By Theorem 2.6, for j > 1, we note that each element <pj in the set is a 
polynomial in A of degree j — 1. Note that we can write any polynomial in A of degree n — 2 
or less as a linear combination of the elements in the set Thus, with respect to 

the spectral measure dp, the quantity <p n is orthogonal to any polynomial with degree n — 2 
or less. Each element p m in the set {p ni } r y~2\ is a polynomial in A of degree m — 1, and 
because m < n we can conclude that p n is orthogonal to <p m with respect to the measure 
dp, i.e. we have 



1 < m < n, 


(4.13) 


where we again suppress the integration limits for notational simplicity. Let us replace dp 
in (4.13) with (dp — dp) + dp and let us replace p n in (4.13) with the right-hand side of 
(2.79). For 1 < m < n we get 


/ / 7°\ 00 / O /o O \ 

(dp - dp) p n Pm + / dp ( <p n Pm J 

n—1 „ n—1 „ 

T ^ ^ I (dp dp) A n j ipj <~p m T ^ ^ / dp (^A n j c pj 


(4.14) 


o 

T 1 m 


Let us define 

Gnm ■■= J (dp - dp) Pn p m . (4.15) 

We will apply the orthonormality given in (4.4) to the second and fourth terms on the 
right-hand side in (4.14) and use (4.15) in the first and third terms on the right-hand side 
in (4.14). This yields 


n—1 n—1 

Gnm T 0 T ^ ^ A n j Gj m T ^ ^ Anj $jm~ 0, 1 T fTl Tl. 

3 =1 3 =1 


(4.16) 


Simplifying (4.16) we obtain the linear system of Gel’fand-Levitan equations 

n — 1 

A nm , + G nrn + ^ A n j Gj m = 0, 1 < Ul < Tl. (4-17) 

3 = 1 


The solution to the inverse problem of recovery of the potential V via the Gel’fand- 
Levitan method is achieved as follows. Given the absolute value of the Jost function /o 
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and the bound-state data, i.e. the zeros of the Jost function f 0 and the corresponding 
Gel’fand-Levitan bound-state norming constants C s appearing in (4.7), with the help of 

(4.1) and (4.6) we form the spectral measures dp and dp. We then evaluate G nm given 
in (4.15). Next, we solve the Gerfand-Levitan system (4.17) and recover A nm . Then, we 
use the values of A n ( n _i) for n — 2, 3,..., b and recover the ordered set {Vi, V 2 , ■ ■ ■, 14} of 
potential values via (2.84). 


Note that the derivation of the discrete Gel’fand-Levitan system (4.17) is similar to 
the following derivation of the continuous version of the Gel’fand-Levitan equation for 

(1.2) with the Dirichlet boundary condition (2.3). Let tp(k,x) be the regular solution to 

(1.2) satisfying (2.15). Let us use !p(k,x) to denote the regular solution to (1.2) with the 
boundary condition (2.15) when V(x) = 0. In fact, we have 


<p( k >x) = 


sin (kx) 


x > 0, 


where we recall that A and k are related to each other as k — \/A. In the continuous case, 
the spectral measures dp and dp are related to the orthonormality relations given by 


dp tp{k,x)tp{k,y) = S(x — y), x, y e [0, + 00 ), 


(4.18) 


J dp [<p(k,x)<p(k,y)] = S(x-y), x, y e [0, + 00 ), (4.19) 

which are the analogs of (4.4) and (4.5), respectively. Note that we use 5(x) to denote the 
Dirac delta distribution. As an analog of (4.15) let us define 


G(x,y) 


(dp-dp) <p(k,x)<p(k,y) 


(4.20) 


We have the analog of the expansion in (2.82), namely 


<p(k, x) = <^(/c, x) + / dz A(x, z) <p(k, z). (4-21) 

Jo 

From (4.19) and (4.21) we conclude that, for y < x, the regular solutions <p(k,x) and 
<p(k,y) are orthogonal with respect to the spectral measure dp, i.e. we have the analog of 
(4.13) given by 


dp <p(k,x)<p(k,y) 


-0, 


y < x. 


(4.22) 


35 



Using (4.21) in (4.22) we obtain the analog of (4.14), namely for y < x we get 

0 = J{dp- dp)y{k,x)y>(k,y) + J clp <p(k,x)<p(k,y) 

+ J dz J(dp — dp) A(x, z) <p(k, z) ip(k,y) + J dz J dp A(x, z) <p(k, z) ip(k, y). 

(4.23) 

Since we use x > y, the second of the four terms on the right-hand side in (4.23) vanishes 
as a result of (4.18). The fourth term there, with the help of (4.18) is seen to be equal 
to A(x,y). Furthermore, using (4.20) in the first and third terms, from (4.23) we get the 
Gel’fand-Levitan equation [5,10,17,18] 


A(x, y) + G(x, y) + / dz A(x, z) G(z, y) = 0, x > y. 


(4.24) 


which is the analog of (4.17). The potential V is recovered from the solution to (4.24) as 

dA(x, x ) 


V(x) = 2 


dx 


dp — \ 


S = 1 

1 y/X 


which is the analog of (2.84). Let us finally remark that the analog of (4.6) in the continuous 
case is [5,10,17,18] 

( * 

AA < 0, 

(4.25) 

. — c/A, A 7 s 0. 

I * |/(Vto)| 2 

where f(V A, 0) with k \/A is the Jost function f(k, 0) and is obtained from the Jost 
solution f(k,x ) appearing in (2.17) by putting x — 0. In (4.25) we assume that there are 
N bound states occurring at A = p, s for s = 1,..., N. The analog of (4.1) in the continuous 
case with the zero potential is then given by [5,10,17,18] 

0, A < 0, 

yx 


dp = 


7T 


c/A, 


A > 0. 


5. TRANSMISSION EIGENVALUES 

In this section we introduce the transmission eigenvalues for the discrete Schrodinger 
equation (1.1) with the Dirichlet boundary condition (2.2). We relate the transmission 
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eigenvalues to the zeros of the key quantity D defined in (5.2) and also equivalently to the 
zeros of the quantity E given in (5.7). We express D in terms of the spectral parameter 
^ and relate it as in (5.8) to the Jost function fo and its relative go, which allows us to 
obtain the relationship (5.10) relating E to fo and go- We introduce the so-called unusual 
case and characterize it in several equivalent ways. We recognize that (5.10) is a discrete 
analog of a Riemann-Hilbert problem, which enables us to solve it uniquely in the usual 
case and recover fo from E. Thus, we uniquely recover the Jost function fo from the 
set of transmission eigenvalues in the usual case. We are then able to use either of the 
Marchenko method and the Gel’fand-Levitan method to uniquely recover the potential in 
class Ab when the transmission eigenvalues are given as input in the usual case. In the 
unusual case, the unique recovery may or may not be possible. This is illustrated with 
some examples in Section 6. 

The transmission eigenvalues [22] for (1.1) with the potential V in class At, and with 
the Dirichlet boundary condition (2.2) are the A-values for which we have nontrivial solu- 

O 

tions ip n and ip n to the system 

V^n+l T 2 1p n 1p n —\ T WiV’n ^Pnt ^ ^ 1> 

o o o o 

~i>n +1 + tyn - V’n-l = A'i/’n, U > 1, 

1 „ (5-1) 

Vo = 0, Vo = o, 

o o 

K ip b = ip b , ipb+l = ipb+l- 

In other words, at a transmission eigenvalue A, we must have (1.1) satisfied for a nontrivial 
wavefunction ip n , the unperturbed problem corresponding to (1.1) with V n = 0 must have 

O O 

a nontrivial solution ip n , the solutions ip n and ip n both must vanish at n = 0, and the 

O 

solutions ip n and tp n must agree at n = b and also at n — b + 1, where b is the positive 
integer appearing in Definition 1.1 and related to the support of the potential V. We note 
that the discrete transmission-eigenvalue problem described in (5.1) is the analog [3] of the 
transmission-eigenvalue problem for the Schrodinger equation with the Dirichlet boundary 
condition (2.3). 

From the first line of (5.1) and the first equation in the third line, it follows that ip n 
must be linearly dependent on the regular solution appearing in (2.42). Similarly, from 
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o 


the second line of (5.1) and the second equation in the third line, it follows that if n must be 
linearly dependent on the regular solution ip n appearing in (2.23). The fourth line of (5.1) 


Wh 

is equivalent to the statement that the column vectors 0 

m+i J 

dependent and hence the matrix formed by using those two vec 


and 


are linearly 


ipb 

V’b+i. 

tors as columns must have 


zero determinant. Consequently, the transmission eigenvalues are exactly the A-values for 


which the two column vectors 

lTb+l J 

the transmission eigenvalues correspond to t 


Tb 


and 


<Pb 

Tb+l 


are linearly dependent. Equivalently, 


determinant via 


D := 


<Pb 


le zeros of the quantity D defined as a matrix 

(5.2) 


<Pb 


I Tb+i Tb+i I 

When V n = 0, from (5.1) it follows that any A is a transmission eigenvalue, and in this case 
from (5.2) we get D = 0. When V n = 0 for n > 2 and V\ ^ 0, from (5.1) it follows that no 
A-value can be a transmission eigenvalue, and in this case from (5.2) we get D = V\. Thus, 
in the analysis of the inverse problem with transmission eigenvalues we can assume that 
b > 2, where b is the positive integer related to the support of V given in Definition 1.1. 


Theorem 5.1 Assume that the potential V belongs to class Ab with 14 ^ 0 and further 
suppose that b > 2, where b denotes the positive integer related to the support of the 
potential. Then, the transmission eigenvalues for (5.1) correspond to the zeros of the 
quantity D defined in (5.2). The quantity D is a polynomial in A of degree 2b — 2, and we 
have 

2b—2 

D=Y, D s^ s , (5.3) 

s =0 

where the coefficients D s for s = 0,1, ..., 2b — 2 are uniquely determined by the potential 
V. In particular, we have 


T>2b-2 — 14, T>2&-3 


b-1 

-V b |4(f>-l) + £Vj 

3 = 1 


(5.4) 


PROOF: The fact that the transmission eigenvalues correspond to the zeros of D has 
already been established in the paragraph containing (5.2). Using (2.70)-(2.72) and (2.76)- 
(2.78) in (5.2), we establish the remaining results stated. | 
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With the help of Theorem 5.1 we obtain the following result. 


Theorem 5.2 Assume that the potential V belongs to class Ab with V& ^ 0 and further 
suppose that b > 2, where b denotes the positive integer related to the support of the 
potential. Then, the quantity D defined in (5.2), as a function of X, has exactly (2b — 2) 
zeros and hence the number of transmission eigenvalues for (5.1) is (2b — 2). In terms of 
the transmission eigenvalues Xi,..., X 2 b- 2 , we have 

2b—2 

B = U n (A(5.5) 

i=i 


2b—2 


Y. x i 


6-1 

4(f>-l) + £Vj. 


1 = 1 


(5.6) 


Hence, the knowledge of the transmission eigenvalues is equivalent to the knowledge of the 


quantity E defined as 


E := 


D 

Vb 


and the transmission eigenvalues correspond to the zeros of E. 


(5.7) 


PROOF: By Theorem 5.1 we know that the zeros of D correspond to the transmission 
eigenvalues for (5.1) and that D is a polynomial in A of degree 2(6 — 1) with the leading 
term 14 A 26-2 . Thus, we have established (5.5). Comparing the coefficients of A 26-3 in (5.3) 
and (5.5), with the help of (5.4) we get (5.6). | 

Recall that we can also use z as the spectral parameter instead of A, where z is related 
to A as in (2.6). Next, we express the key quantity D given in (5.2) in terms of the spectral 
parameter 2 . This will help us to solve the relevant inverse problem to recover the potential 
with the help of the Marchenko method or the Gel’fand-Levitan method. 

Theorem 5.3 Assume that the potential V belongs to class Ab with 14 ^ 0 and further 
suppose that b > 2, where b denotes the positive integer related to the support of the 
potential. Then: 

(a) The key quantity D given in (5.2) can be expressed in terms of the spectral parameter 
z via 

D = ^ _ 1 ^_ 1 (fo - go ), (5.8) 
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where fo is the Jost function appearing in (2-46) and go is the quantity in (2-49). The 
expression (5.8) holds for h — 1 as well. 

(b) The quantity (fo — 1 )/V b is uniquely determined by the transmission eigenvalues for 
(5.1). The quantity Koi/Vb is also uniquely determined by the transmission eigenval¬ 
ues, where K 0 \ Is the coefficient of z in the polynomial expansion (2-46) in terms of 
z for the Jost function fo- In fact, the value of Koi/V b is equal to the coefficient of z 
when (fo — 1 )/V b is expressed in z as a polynomial. 


(c) 


The identity (5.6) is equivalent to 

26-2 

E = - 1 ) + 

3 = 1 



(5.9) 


(d) Unless the quantity Koi/V b is equal to 1, where the case Koi/Vi, = 1 corresponds to 
the unusual case, V), is uniquely determined by the transmission eigenvalues by solving 
(5.9) for V b . 

(e) Unless we are in the unusual case, the transmission eigenvalues uniquely determine the 
Jost function fo- Thus, unless we are in the unusual case, the transmission eigenvalues 
uniquely determine the ordered set {V \,..., V b ) of potential values. 

PROOF: Using (2.21) and (2.42) in the first and second columns, respectively, of (5.2) we 
establish (a). Using (5.7) in (5.8), we get 

(5.10) 

From Theorem 5.2 we know that the left-hand side of (5.10) is uniquely determined by the 
transmission eigenvalues of (5.1). From Theorem 2.2 we know that fo — 1 is a polynomial 
in z of degree (2b — 1) and the quantity <7o — 1 is a, polynomial in z~ x of degree (2b — 1). 
Thus, given the left-hand side of (5.10), we can uniquely determine (fo — 1)/V& as follows. 
From Theorem 2.2 it follows that (/ 0 — 1)/V& is a “plus” function and (go — 1 )/V b is a 
“minus” function, and hence writing the left-hand side of (5.10) in the spectral parameter 
z, we see that the terms containing the positive powers of z make up (fo — 1 )/H and the 
terms containing the negative powers of z make up (go — 1) /V b . Having (fo — 1) /V b at hand, 
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we can use (3.17) with n — 0 and m = 1 to obtain the value of K m /Vi,. We remark that 
Koi/Vb is the coefficient of z when (fo — 1 )/V\ ) is written in terms of z as a polynomial in 
z. Thus, we have proved (b). From the second equality in (2.47) we get 


K , 


01 


1 


6-1 


V b 


1 + H^> 

1 = 1 


(5.11) 


and hence using (5.11) in (5.6) we obtain (5.9), which establishes (c). By (b) we know that 
K 01 /V b is uniquely determined by the transmission eigenvalues and hence we can solve 
(5.9) for Vb and determine 14 uniquely in terms of the transmission eigenvalues provided 
Koi/Vb ^ 1. Thus, we have also established (d). From (b) we already know that the 
transmission eigenvalues uniquely determine (fo — 1)/Vj,, and from (d) we know that in the 
usual case the transmission eigenvalues uniquely determine V b . Thus, in the usual case the 
transmission eigenvalues uniquely determine fo- By Theorem 3.1 we know that fo uniquely 
determines the potential V. Thus, the proof of (e) is complete. | 


The following result is needed to prove Theorem 5.5. It shows that a transmission 
eigenvalue for (5.1) cannot occur at a zero of the Jost function / 0 . 

Proposition 5.4 Assume that the potential V belongs to class Ab specified in Defini¬ 
tion 1.1. Let fo be the Jost function appearing in (2.f6), go the quantity in (2.f9), and D 
the quantity defined in (5.2). Let z be the spectral parameter appearing in (2.12). Then, 
the quantity D and fo cannot vanish simultaneously at any z-value. 


PROOF: From (5.8) we see that if D and fo vanished at the same 2 -value with z ±1, then 
we would also have go vanishing at that z- value, which would contradict Theorem 2.4(e). 
Thus, D and fo cannot simultaneously vanish at any z- value with z LI. On the other 
hand, using (2.59) and (2.61) with n — 0 as well as (2.41) and (2.64) in (5.8), we get the 
expansion for D around z = 1 as 

D(z)= fo(l) + 0(z-l), z^ linC, (5.12) 

where we recall that / 0 (1) denotes the derivative of fo with respect to z evaluated at z — 1. 
Because of (2.65) we know that / 0 (1) ^ 0 if /o(l) = 0. Thus, D and fo cannot simulta¬ 
neously vanish at z — 1. A similar argument shows that D and fo cannot simultaneously 
vanish at z = — 1 either. | 
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From (2.20) and (2.45) it follows that the scattering matrix for the unperturbed system 
is given by 

O 

5=1. 

The following theorem explores the connection between the transmission eigenvalues and 
the energies at which the scattering from the perturbed system (1.1) and from the unper¬ 
turbed system (2.5) coincide. It shows that at a transmission eigenvalue the scattering 
matrix S of the perturbed system takes the value of 1. With the possible exception of 
A = 0 and A = 4, it shows that any A-value at which S becomes equal to 1 is a transmis¬ 
sion eigenvalue. It also gives some necessary and sufficient conditions for each of A = 0 
and A = 4 to be a transmission eigenvalue for the system (5.1). 

Theorem 5.5 Assume that the potential V belongs to class Ab with 14 ^ 0 and further 
suppose that b > 2, where b denotes the positive integer related to the support of the 
potential. Then: 

(a) At a transmission eigenvalue A for (5.1), the scattering matrix S given in (2.f5) takes 
the value 1. Thus, at a transmission eigenvalue A the scattering from the perturbed 
system (1.1) with the Dirichlet boundary condition (2.2) agrees with the scattering 
from the unperturbed system (2.5) with the Dirichlet boundary condition (2.2). 

(b) With the possible exception of z — 1 and z — —1, corresponding to X — 0 and A = 4, 
respectively, each X-value at which the scattering matrix S takes the value of 1 is a 
transmission eigenvalue of (5.1). 

(c) If the scattering matrix S takes the value 1 at X — 0, then A = 0 is a transmission 
eigenvalue if and only if fo is nonzero and fo is zero at X — 0. 

(d) If the scattering matrix S takes the value 1 at X — 4, then X — A is a transmission 
eigenvalue if and only if fo is nonzero and fo is zero at A = 4. 

(e) If the scattering matrix S takes the value 1 at X — 0, then A = 0 is a transmission 
eigenvalue if and only if S — 1 has a double zero at A = 0. 

(f) If the scattering matrix S takes the value 1 at X — 4, then A = 4 is a transmission 
eigenvalue if and only if S — 1 has a double zero at X — A. 
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PROOF: Recall that A = 0, as a result of (2.6), corresponds to z — 1 in a one-to-one 
manner and that A = 4 corresponds to z = — 1 in a one-to-one manner. With the help of 
(2.45), let us write (5.8) as 

iz-z~ 1 )^- = l-S. (5.13) 

Jo 

By Theorem 5.1 we know that the transmission eigenvalues correspond to the zeros of D. 
First, consider a zero of D occurring at a 2 -value with z ^ ±1. From Proposition 5.4 we 
know that /o must be nonzero at such a z- value. Hence, the left-hand side of (5.13) is zero 
at that value, causing S to take the value of 1 at that 2 -value. If the zero of D occurs at 
2 = 1, then from Proposition 5.4 we know that fo\ z =i 0 and hence the left-hand side of 
(5.13) is zero, which results in S taking the value of 1 at 2 = 1. A similar argument shows 
that if 2 = —1 corresponds to a zero of D, then we have S taking the value of 1 at 2 = —1. 
Thus, the proof of (a) is complete. The proof of (b) is given as follows. If S = 1 at some 
2 -value other than 2 = ±1, then the left-hand side of (5.13) must be equal to zero at that 
2 -value, and hence (5.13) implies that D must vanish at that same 2 -value. Thus, that 
2 -value, being a zero of D must be a transmission eigenvalue, establishing (b). To prove 
(c) we proceed as follows. Using (2.59) and (2.61) in (2.45), we can expand the scattering 
matrix S around 2 = 1 as 

foil) - (z - 1) foil) + 0(iz- l) 2 ) 


/„(l) + ( 2 -l)/ 0 (l) + O((2-l) 2 )’ 
which yields, for 2 -> 1 in C, 

7o(i) 


2 —» 1 in C, 


1 — 2(2 — IITItt + C ((2 — l) 2 ), / 0 (1) # 0, 

S(z) = { Ml) 

-1 + 0 ( 2 - 1 ), /»(!) = 0 . 


(5.14) 


Thus, we have F(l) = 1 if / 0 (1) ^ 0 and we get F(l) = —1 if /o(l) = 0. On the other hand, 
from (5.12) we know that D vanishes at 2 = 1 if and only if /o(l) = 0. Therefore, at 2 = 1 
we have S = 1 and D — 0 if and only if /o(l) ^ 0 and /o(l) = 0. Hence, the proof of (c) is 
complete. The proof of (d) is obtained in a similar way as in the proof of (c). From (5.14) 
we see that S — 1 has a double zero at 2 = 1 if and only if we have foj^O and Jo = 0 at 
2 = 1. Thus, (e) is equivalent to (c). Similarly, one can show that (f) is equivalent to (d), 
completing the proof of our theorem. | 
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Let us now investigate the unusual case in the transmission-eigenvalue problem for 
(5.1). Recall from Theorem 5.3(d) that the unusual case occurs for (5.1) when 

26—2 

Y A,- = 4(6 - 1), (5.15) 

3 = 1 

where 6 is the integer related to the support of the potential V and the A 7 -values for 
j — 1,..., 26—2 denote the transmission eigenvalues. From (5.9) and (5.11), or equivalently 
from (5.6) and (5.15), we see that the unusual case can also be identified as the case with 

6-1 

Y V i=°- < 8 ' 16 > 

3 = 1 

In the next theorem we characterize the unusual case in several different ways. 

Theorem 5.6 Assume that the potential V belongs to class Ab with 14 ^ 0 and further 
suppose that b > 2, where b denotes the positive integer related to the support of the 
potential. Then, the unusual case for the transmission eigenvalue problem for (5.1) occurs 
if and only if any one of the following equivalent conditions is satisfied: 

(a) The transmission eigenvalues A j satisfies (5.15). 

(b) The potential V satisfies (5.16). 

(c) The coefficient K (n of z in the polynomial expression in z given in (2-46) for the Jost 
function fo is equal to 14 - 

(d) The coefficient Ko( 2 b- 2 ) of z 2b ~ 2 in the polynomial expression in z given in (2-46) for 
the Jost function fo is zero. 

(e) The coefficients of z and z 2b ~ x in the polynomial expansion in z for (fo — 1)/14 ore 
both equal to 1. 

(f) The coefficient of z 2b ~ 2 in the polynomial expansion in z for (fo — 1)/14 Is zero. 

PROOF: Note that (a) can be used as the definition of the unusual case. From (5.9), (5.15), 
and (5.16) we get the equivalence of (a) and (b). Using (5.16) in the second equality in 
(2.47) we see that the unusual case occurs if and only if Kqi = 14, establishing (c). Using 
(5.16) in the second equation in (2.48) we observe that the unusual case is equivalent 
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to having Ko( 2 b- 2 ) — 0; which establishes (d). From (2.46)-(2.48) it follows that the 
coefficients of z and z 26-1 in the polynomial expansion in z for (/o — 1)/V& are given by 
Kqi/Vi, and iF 0 ( 26 -i)/H) respectively, and hence with the help of (5.11) and (5.16) we see 
that Kqi/Vi, — 1 in the unusual case and K G r 2 b-i)/Vb = 1 always, which establishes (e). 
From (2.46) and the first equality in (2.48), or equivalently from (2.50), we see that the 
coefficient of z 2b ~ 2 in the polynomial expansion in z for (/q — 1 )/V, is zero if and only if 
(5.16) holds, establishing (f). | 

In the usual case we already know that the transmission eigenvalues uniquely de¬ 
termine the potential. By presenting some explicit examples, we will show in the next 
section that the transmission eigenvalues for (1.1) may or may not uniquely determine the 
potential V in the unusual case. We will show that in the unusual case we may have a 
one-parameter family of real-valued potentials corresponding to a given set of transmission 
eigenvalues, we may have a finite number of distinct real-valued potentials corresponding 
to a given set of transmission eigenvalues, or we may have a unique real-valued potential 
corresponding to a given set of transmission eigenvalues. 

6. EXAMPLES 

In this section we present various explicit examples to illustrate the results of the 
previous sections. In particular, we illustrate the recovery of the potential from the trans¬ 
mission eigenvalues, the Marchenko method, the Gel’fand-Levitan method, the unusual 
case, and the bound states. 

In the first example below we illustrate the fact that the bound-state information is 
contained in the scattering matrix S when the corresponding potential V belongs to class 
Ab described in Definition 1.1. As stated in Theorem 2.5(e) the bound-state energies cor¬ 
respond to the poles of the scattering matrix, and as indicated in (3.19) the corresponding 
Marchenko bound-state norming constants are obtained from S via a residue evaluation. 

Example 6.1 Assume that the potential is supported at n — 1 with V\ ^ 0, and hence 
V n — 0 for n > 2. From (2.24) we have f n — z n for n > 1. Using (2.29) and (2.46) we 
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obtain the Jost function as 


fo = V\z + 1. 


( 6 . 1 ) 


Thus, the only zero of fo occurs at z — —1/V\. From Theorem 2.5(a) we know that a bound 
state exists if the zero of fo is in one of the real line segments z G (0,1) and z G (—1,0). 
Thus, there is exactly one bound state if V\ < —1 or V\ > 1. Assume that one of these 
conditions is satisfied so that we have a bound state at z — —1/V\. From (6.1) we then get 

/.L-i/vi = (-£)"• 

Using (6.2) in (3.15) we evaluate the Marchenko bound-state norming constant as 


Cl 


1 



which simplifies to 

ci “A?- 1 - 

From (6.1), replacing z by z~ 1 , we get the quantity go appearing in (2.49) as 


(6.3) 


go — Ur 1 +1- 


(6.4) 


Using (6.1) and (6.4) in (2.45), we obtain the scattering matrix S as 

c _ V lZ ~ l + 1 
V lZ + l 


(6.5) 


Let us verify that (3.16) yields the Marchenko bound-state norming constant given in (6.3). 
From (6.5) we evaluate the residue of S/z at z = — 1/Ui, and we have 


Res 


~S _J_' 
7’ Vi_ 



and a comparison with (6.3) shows that (3.16) is verified. Regarding the transmission 
eigenvalues, from (5.3) we get D = \\ and hence there are no transmission eigenvalues. 


In the second example given below, we demonstrate that the transmission eigenvalues 
may have multiplicities greater than one. We also illustrate that in the unusual case there 
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may be infinitely many distinct potentials corresponding to a given set of transmission 
eigenvalues. 

Example 6.2 Let us assume that V n = 0 for n > 1 except when n — b for some fixed 
positive integer b and that 14 is a nonzero real parameter. By Theorem 5.6(b) we see 
that this corresponds to a special case of the unusual case. Proceeding as in the proof of 
Theorem 2.6 we can determine all the values of the regular solution tp n . We obtain 


Vb ~ Vbi Vb+l — Vb +1 + Vb ipb, 


( 6 . 6 ) 


where we recall that ip n is the regular solution when V n = 0 and it appears in (2.21) and 
(2.22). Using (6.6) in (5.2), we simplify the determinant defining D and get 


D = 


ifib 0 

° TZ ° 

< Pb+1 Vb ^b 

From (5.7) and (6.7) we observe that the quantity E is given by 

2 


Using (2.22) in (6.8) we have 


(6.7) 


E = [Vb 


E = 


sin 2 jbO) 
sin 2 6 


( 6 . 8 ) 


(6.9) 


Comparing (5.5) and (6.9) we see that there are 2(6 — 1) transmission eigenvalues A j 
including multiplicities, and in fact the transmission eigenvalues occur when 6 = jit/b for 
j = 1, 2,..., 6 — 1 and each has multiplicity two. Using (2.11) we obtain the A-values for 
the transmission eigenvalues as 


A 7 ’ — 2 — 2 cos I 

' 6 


J = 1, 2,.. •, 6 — 1. 


( 6 . 10 ) 


Thus, with the help of (5.5), (5.7), and (6.10) we can write (6.9) as 

6-1 


£ =n 

j =i L 


A — 2 + 2 cos 


J7T 


( 6 . 11 ) 


From (6.11) we observe that there are no transmission eigenvalues if 6 = 1, there is a trans¬ 
mission eigenvalue at A = 2 with multiplicity two if 6=2. and there are the transmission 
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eigenvalues at A = 1 and A = 3 with multiplicities two if b — 3. It is clear from (6.11) that 
the value of 14 cannot be determined from the transmission eigenvalues, and hence there 
is a one-parameter family of potentials corresponding to the transmission eigenvalues for 
b> 2, with 14 being the nonzero real parameter. 

In the third example given below we illustrate that the transmission eigenvalues may 
be all real, may be all complex, or may have multiplicities greater than one. 

Example 6.3 Assume that b — 2 in Definition 1.1 and hence V n = 0 for n > 3. Further 
assume that 14 4 0. Thus, the potential is known by specifying 14 and 14- From (2.24) 
and (2.29)-(2.31) we get f n — z n for n > 2 and 

fi = V 2 z 2 + z, f 0 = 14 z 3 + VlV 2 z 2 + (14 + 14 )z+ 1. (6.12) 

Note that we can get / 0 given in (6.12) also from (2.50). Thus, the quantity g 0 appearing 
in (2.49) is obtained from /o given in (6.12) by replacing z with z -1 . We have 

go = 14 W 3 + ViV 2 z~ 2 + (14 + 14) W 1 + 1. (6.13) 

Using (6.12) and (6.13) in (5.8) we obtain 

n _ 14 O? 3 - £- 3 ) +141404 - Z~ 2 ) + (14 +14) (z - z- 1 ) 

z-z - 1 

or equivalently 


D = V 2 (z 2 + 1 + z~ 2 ) + V^iz + z- 1 ) + (14 + 14). 
Using (2.8) in (6.14) we express D as a function of A as 

D = 14A 2 - 14(4 + 14 )A + (14 + 21414 + 414). 
Thus, the quantity E is obtained by using (6.15) in (5.7) as 

E = A 2 - (4 + 14)A + ^4 + 214 + . 


(6.14) 


(6.15) 


(6.16) 
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We know by Theorem 5.2 that the transmission eigenvalues correspond to the zeros of E 
given in (6.16). Thus, from (6.16) we see that the two transmission eigenvalues are given 
by 


A 



(6.17) 


Hence, depending on the potential, we may have two real distinct transmission eigenvalues, 
a double real transmission eigenvalue, or a pair of complex conjugate transmission eigenval¬ 
ues. For example, we have a double transmission eigenvalue at A = 0 if (Vi, V 2 ) = (—4, —1), 
we have (Ai,A 2 ) = (i, — i) when (Ui, V 2 ) = (—4,—4/5), we have (Ai,A 2 ) = (3,-3) when 
(Vi, V 2 ) = (—4, 4/5), and we have (Ai, A 2 ) = (1 + i, 1 — i) when (Vf, V 2 ) = (— 2 , — 1 ). From 
Theorem 5.6(b) it follows that the unusual case occurs when V\ = 0, in which case we have 


E = (A — 2 ) 2 . 


As already explained in Example 6 . 2 , in the unusual case, corresponding to the transmission 
eigenvalues (Ai, A 2 ) = ( 2 , 2 ) we have a one-parameter family of potentials with V\ = 0 and 
V 2 being the nonzero real parameter. 

In the fourth example given below, we illustrate Theorem 5.5(c) and Theorem 5.5(d). 

Example 6.4 As in Example 6.3, let us assume that b — 2 with V 2 7 ^ 0, and hence V n — 0 
for n > 3. Let us first illustrate Theorem 5.5(c). Using / 0 given in (6.12), at z — 1 we 
obtain 

/o(l) = 2 U 2 + Ui + WU 2 + 1 , /o(l) - 4U 2 + Vi + 2 ViV 2 . (6.18) 

From (6.18) we see that /o(l) — 0 is equivalent to 


V 2 = 


Vi 

4 + 2Vi ’ 


(6.19) 


in which case from (5.14) we get 5'| z =i = 1. Using (6.19) in the first equality in (6.18), we 
see that /o(l) 7 ^ 0 is equivalent to V\ 7 ^ —2. From the denominator on the right-hand side 
in (6.19) we already see that we must have Vf 7 ^ —2 for (6.19) to make sense. Let us now 
determine when A = 0 becomes a transmission eigenvalue. From (6.14) we see that A = 0 
is a transmission eigenvalue if and only if the sum of the last three terms on the right-hand 
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side in (6.14) is zero, i.e. 


4 + 2Vi + -^ = 0, 

V2 

which is equivalent to having (6.19). Thus, A = 0 is a transmission eigenvalue if and 
only if we have /o(l) ^ 0 and /o(l) = 0, verifying Theorem 5.5(c). Let us now illustrate 
Theorem 5.5(d). Evaluating /o of (6.12) at z — —1 we get 

/o(-l) = -2E 2 -1A + ViV 2 + 1, /„(-!) = 4V 2 + V 1 - 2V 1 V 2 . (6.20) 


From (6.20) we see that /o(—1) = 0 is equivalent to 


U 2 = 


Vi 

4 - 2 VI' 


( 6 . 21 ) 


Using (6.21) in the first equality in (6.20), we see that /o(—1) ^ 0 is equivalent to Vi ^ 2. 
From the denominator on the right-hand side in (6.21) we already see that we must have 
Vi ^ 2 for (6.21) to make sense. Let us now determine when A = 4 becomes a transmission 
eigenvalue. From (6.17) we see that A = 4 is a transmission eigenvalue if and only if 


4 




or equivalently if and only if 




YL_Yjl 

4 V 2 ' 


( 6 . 22 ) 


Simplifying (6.22) we see that (6.22) is equivalent to (6.21). Thus, A = 4 is a transmission 
eigenvalue if and only if we have /o( —1) ^ 0 and /o(—1) = 0, verifying Theorem 5.5(d). 


In the fifth example given below we illustrate that, in the unusual case, we may have 
a unique real-valued potential V in class At, corresponding to a given set of transmission 
eigenvalues, we may have a finite number of distinct real-valued potentials corresponding to 
a given set of transmission eigenvalues, or we may have infinitely many distinct real-valued 
potentials corresponding to a given set of transmission eigenvalues. 


Example 6.5 Assume that b — 3 with V 3 7 ^ 0, and hence V n = 0 for n > 4. By Theorem 5.2, 
the system (5.1) has four transmission eigenvalues uniquely determined by the ordered set 
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{Vi,V 2 ,Vs}. By Theorem 5.6(b) we know that the unusual case occurs when V> — —V\. 
Let us see to what extent the transmission eigenvalues determine the potential V in the 
unusual case. Proceeding as in the proof of Theorem 2.6, we can evaluate the values of 


the regular solution ip n . We get 

<p 3 = A 2 - 4A + 3 - Vi 2 , (6.23) 

<P4 = -A 3 + (6 + u 3 )A 2 + (Vi 2 - 4V 3 - 10)A + (4 - Vi - 2Vf + 3E 3 - V 3 V?), (6.24) 

and hence we also have the values of the regular solution ip n corresponding to V n = 0, and 
in particular from (2.75) we obtain 

<p 3 = A 2 -4A + 3, <p 4 = -A 3 + 6A 2 - 10A + 4. (6.25) 

Since b — 3 in our example, using (6.23)-(6.25) in (5.2), with the help of (5.7) we get 

E = A 4 - 8A 3 + (22 - 7 )A 2 + (4 7 + e - 24)A + (9 - 3 7 - 2e), (6.26) 


where we have let 


7 := W + 


Vi 



(6.27) 


Recall that we assume Vi = —Vi, which signifies the unusual case. Since the inverse 
problem of determining the potential V from the transmission eigenvalues is equivalent to 
the recovery of V from E, our inverse problem is equivalent to the determination of Vi 
and V 3 from the ordered set { 7 , e}. Let us analyze this problem algebraically and see to 
what extent we can determine the real parameter Vi and the nonzero real parameter V 3 
from the ordered set { 7 , e}. Eliminating V 3 in (6.27) we obtain 


Vf- 7 VA + e = 0. 


(6.28) 


Given the ordered set { 7 , e}, the solution to the algebraic equation in (6.28) may yield 
one, two, or three distinct real values for Vi. Once we obtain Vi by solving (6.28), we can 
recover V 3 from the second identity in (6.27). If 7 = 0 and e ~ 0. then Vi = 0 is the 
solution with multiplicity three. This special case is analyzed in Example 6.2 already, and 
it corresponds to the unusual case. We already know that in this unusual case, V 3 cannot 
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be determined from the knowledge of transmission eigenvalues. In this particular case the 
four transmission eigenvalues are given by 

Ai = 1, A 2 = 1, A 3 = 3, A 4 = 3, 

corresponding to the one-parameter family of potentials with V n — 0 for n ^ 3 and V 3 
being any nonzero real constant. Let us next consider the case 7 = 0 and e 7 ^ 0. In this 
case from (6.27) and (6.28) we get 

Vi = -tfi, v 3 = ^, 

indicating that the transmission eigenvalues uniquely determine the potential V in this 
particular case. For example, the choice 7 = 0 and e = 1 yields the unique potential V 
with 

(Vi, v 2 , V 3 ) = (-1,1,1), 

corresponding to the transmission eigenvalues 

Ai = 1.47511, A 2 = 0.509784, A 3 = 3.00755 + 0.513116i, A 4 - 3.00755 - 0.513116 i, 

which are obtained as the zeros of E given in (6.26). Note that we use an overline on a 
digit to indicate a round off. Let us next consider the case where (6.28) has three distinct 
real solutions for Vi. For example, when 7 = 7 and e = 6, we have three distinct solutions 
for V\ in (6.28), namely, 1, 2, and —3. In this case, the transmission eigenvalues A j are 
obtained as the zeros of E from (6.26) as 

Xi — 4, A 2 = 3.85577, A 3 = 1.32164, A 4 = -1.17741, 

corresponding to exactly the three distinct potentials with the values of (Vf, V 2 . V 3 ) given 
by 



When 7 = 3 and e = 2, we get exactly two distinct real V\ -values as solutions to (6.28), 
namely Vi = 1 and V\ = —2. In this particular case, the transmission eigenvalues A ? are 
obtained as the zeros of E from (6.26) as 

Ai = 1.2495, A 2 = -0.258627, A 3 = 3.50456 + 0.309892 i, A 4 - 3.50456 - 0.309892 i, 
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corresponding to exactly the two distinct potentials with the values for (Vi, V 2 , V 3 ) given 
by 

(l,-l,|), (- 2 , 2 , 2 ). 

As a summary, in this example, corresponding to the four transmission eigenvalues, we have 
demonstrated that we may have a unique potential in class At,- two distinct potentials, 
three distinct potentials, or a one-parameter family of infinitely many potentials. 


In the sixth example below we illustrate the use of the Marchenko method to obtain 
the potential from the transmission eigenvalues. 


Example 6.6 Assume that 6 = 2 with V 2 ^ 0, and hence V n — 0 for n > 3. We have 
2(6 — 1 ) equal to 2 , and hence there are two transmission eigenvalues, as indicated in 
Theorem 5.2. Let us assume that they are given by 


Ai 


11 + ^57 

4 


A 2 


11 - \/57 
4 


(6.29) 


Using (5.5), (5.7), and (6.29) we see that the quantity E defined in (5.7) is given by 


E= A 


11 + ^ 57 ' 

If - 


A 


11 - ^57 N 

If - 


(6.30) 


Using (2.8) in (6.30) we express E in terms of the spectral parameter 2 and obtain 


2 3 3 , 2 

E — z 2 +-z — l + - z~ l + W 2 . 

Z Z 


Thus, we have 

(z — z~ x ) E — z 3 + ^ z 2 — 2z + 2z~ x — ^ z~ 2 + W 3 . (6.31) 

Z Z 


From (6.31), with the help of (5.10) we extract (/o — 1 )/V 2 by taking the positive powers 
of z and get 


/o-l 

U 2 


= z 6 + -z z - 2z. 
z 


(6.32) 


As seen from Theorem 5.3(b), the coefficient of z on the right-hand side in (6.32) gives us 
K 01 /V 2 , and hence 


K01 

V 2 


(6.33) 


53 



Using (6.29) and (6.33) in (5.9), we write (5.9) as 


- = 4 + (-2 - 1) V 2 , 


and hence recover V 2 as 


V 2 = -- 


(6.34) 


Note that the first equality in (2.48) indicates that the coefficient of z 2 in (6.32) must be 
Vi, i.e. we have 

However, if b were a large integer, we would not be able to obtain the ordered set 
{V,. Vh,.... Vi,- 1 } of potential values in such a straightforward manner. Thus, we will 
use the Marchenko procedure to illustrate the recovery of V\ because it is a more system¬ 
atic way of recovering the potential. Using (6.34) in (6.32) we obtain fo as 


1 3 3 2 

= ~2 Z 4 Z +Z + L 


(6.35) 


From (2.46) and (2.49) we see that the quantity go is obtained by replacing z with z 1 in 
(6.35) and we have 

9o = ~\ z ~ 3 ~ \ z ~ 2 + z_1 + L (6-36) 


Using (6.35) and (6.36) in (2.45) we obtain the scattering matrix S as 

1 -3 3 -2 -l 1 

o “2* *4* + * +1 

1 3 3 2 ' 

~2 Z 4 Z +Z + 1 


(6.37) 


Let us illustrate the determination of the potential V from S given in (6.37). We determine 
the zeros of fo inside the unit circle \z\ — 1, which is equivalent to finding the poles of S, as 
indicated in Theorem 2.5(e). Note that fo has three zeros at z — z a for s = 1,2, 3, where 

1 — a/ 17 l + v^ 

*1 = - 7 -> - 7 -> ^3 = - 2 . 


We have 


Z! = -0.780776, Z 2 = 1.28078, = -2, 
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and hence only z\ corresponds to a bound state based on the criteria given in Theo¬ 
rem 2.5(a), i.e. among the three ^-values, only z\ lies in the union of the intervals (—1, 0) 
and (0, 1). Next, for n > 1 we determine M n defined in (3.14). The integral in (3.14) can 

be evaluated in terms of the residues of S at the two poles z — 0 and z = Z\ inside the 

unit circle T. From (3.14) and (3.16) we see that the summation term in (3.14) consists of 
one term only and it cancels the contribution from the integral due to the pole at z = Z\. 
Thus, the only contribution to M n comes from the residues at z = 0 and we have 

Mi = —Res [S, 0], M 2 — —Res [ zS , 0], M 3 = —Res [ z 2 S , 0] . (6.38) 

Using (6.37) in (6.38) we get 

Mi = — ^ M 2 = i M 3 = i, ( 6 . 39 ) 

M n = 0, n > 4. (6.40) 

Using (6.39) and (6.40) let us write the Marchenko system (3.12) by retaining only the 
nonzero terms. We get 

A01 + Mi + A01M2 + A02M3 = 0, 

A 02 + M 2 + A.01M3 = 0, 

< (6.41) 

A 03 + Ms = 0 , 

„ A 12 + Ms — 0 . 

Thus, from the Marchenko system (3.12) we get K nm = 0 for 0 < n < m, with the 
exception of the four nonzero terms Kqi, A 02 , A' ( ,s, and K\ 2 - By solving the linear system 
(6.41) we obtain 

3 11 

A"oi = 1, A 02 = --, A 03 =: —, A 12 — — • 

Finally, we recover the potential via (2.51) and obtain 

Ul - A01 - A12, U 2 = K 12 - A 23 , 


which yields 



U 2 = 


1 

2 
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In the seventh example below we illustrate the Gel’fand-Levitan method presented in 
Section 4. 

Example 6.7 In this example we illustrate the use of the Gel’fand-Levitan method to 
recover the potential from the transmission eigenvalues. Assume that 6 = 2 and U 2 7 ^ 0 . 
Hence, V n = 0 for n > 3. We have 2(6— 1) equal to 2, and hence there are two transmission 
eigenvalues as asserted in Theorem 5.2. Let us assume that they are given by 


Ai = - 1 , A 2 = 4. (6.42) 

From (5.5) and (6.42) we see that the quantity E defined in (5.7) is given by 

E = A 2 — 3A — 4. (6.43) 

Using (2.8) in (6.43) we express E in terms of the spectral parameter 2 and obtain 

(z - z- 1 ) E = z 3 - z 2 - 5z + 5 W 1 + W 2 - W 3 . (6.44) 


From (6.44), with the help of (5.10) we extract (/o — 1 )/V 2 by taking only the positive 
powers of z and get 

= z 3 — z 2 — 5 z. (6.45) 

V2 


As seen from Theorem 5.3(b), the coefficient of z on the right-hand side in (6.45) gives us 
K 01 /V 2 , and hence we have 


Kqi = 

U 2 


(6.46) 


Using (6.42) and (6.46) in (5.9), we write (5.9) as 


3 = 4 +(-5-1) V 2 , 


and hence recover U 2 as 

U 2 = 1 (6.47) 

6 

Note that the first equality in (2.48) indicates that the coefficient of z 2 in (6.45) must be 
Vi, i.e. we have 

Ui = -1. 
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However, if b were a large integer, we would not be able to obtain the ordered set 
{Vi, V 2 ,..., Vb-i} of potential values in such a straightforward manner. Instead, we will 
use the Gel’fand-Levitan procedure to illustrate the recovery of V\ because that is a fun¬ 
damental procedure to obtain the potential. Using (6.47) in (6.45) we obtain / 0 as 

/„ = U 3 -U 2 + U + 1 . (6.48) 

6 6 6 

From (2.46) and (2.49) we see that the quantity go is obtained by replacing z with z^ 1 in 
(6.48) and hence we have 

9o = U" 3 -U “ 2 + U“ 1 +1. (6.49) 

6 b b 


We determine the zeros of fo inside the unit circle \z\ — 1. Note that /q has three zeros at 
z = z s for s = 1,2, 3, where 


Zl 


-1 - 
2 


-1 + Vl3 

2 


z 3 = 2 . 


(6.50) 


We have 


Z! = -2.30278, = 1.30278, £3 - 2, 


and hence among the three -values none lie in the union of the real intervals (— 1 , 0 ) and 
(0,1). Thus, by Theorem 2.5(a) there are no bound states. Let us now evaluate G nm given 
in (4.15). Since there are no bound states, from (4.3) and (4.12) we obtain 

d(p-°p) = ^-(z-z- 1 f (1--^-)-, ze T+, (6.61) 

2™ v JO , 9 o J z 

which is extended in a straightforward manner from T + to the whole circle T. Using (2.21) 
and (6.51) in (4.15) we get 



where the factor 1/2 in front of the integral in (6.52) is due to the fact that we now integrate 
along the entire unit circle T in the positive direction instead of integrating along the upper 
semicircle T + . We remark the symmetry G nm — G , rvn , as seen from (6.52). The integral 
in (6.52) can be evaluated by using the residues at the poles of the integrand inside the 
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unit circle. Such poles occur at z = 0 and z — 1 /Zj for j = 1, 2, 3, where the ^-values are 
given in (6.50). We only need to determine V\ and V 2 via the Gel’fand-Levitan method, 
and hence it is sufficient to determine G nm only for some small values of n and m. Since 
we will use (2.84) for n — 1 and n — 2, from (2.84) we see that it is enough for us to set 
up the Gel’fand-Levitan system only for 1 < m < n < 3. Using (6.48) and (6.49) in (6.52) 
we obtain 

Gn — 0, G 21 = 1, G 22 — 1, G 31 = 1, G 32 = —. (6.53) 

6 

We write the Gel’fand-Levitan system (4.17) for 1 < m < n < 3 as 

{ ^21 + G 21 + A 21 Gn — 0 , 

^31 + G 31 + A 31 G 11 + A 32 G 21 = 0, (6.54) 

^-32 + G 32 + ^4.31^12 + A 32 G 22 — 0. 

Using (6.53) in (6.54) we get 

A 2 i + 1 = 0 , 

< ^31 + 1 + A 32 = 0 , 

„ 11 

^-32 H p h -4.31 + 4.32 = 0, 

6 

which is uniquely solvable and yielding 

4.21 = —1, 4 3 i = —A 32 = —(6.55) 

6 6 

We recall that 4io = 0 as a result of (2.81). From (2.84) we have 

Ul = 4 2 i — 4i 0 , V 2 = 4 32 — 4 2 i, (6.56) 

and hence using (6.55) and (6.56) we get 

W = -l, V 2 = ^. 

6 
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